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Abstract 

An  electrically-regenerated  electrosorption  process  known  as  carbon  aerogel  CDI  has  been 
developed  by  Lawrence  Livermore  National  Laboratory  (LLNL)  for  continuously  removing  ionic 
impurities  from  aqueous  streams.  A  salt  solution  flows  in  an  unobstructed  channel  formed  by 
numerous  pairs  of  parallel  carbon  aerogel  electrodes.  Each  electrode  has  a  very  high  BET 
surface  area  (2. 0-5. 4x1 06  ft2  lb"1  or  400-1100  m2  g'1)  and  very  low  electrical  resistivity  (<  40  mQ 
cm).  BET  surface  areas  of  1.3xl07  ft2  lb'1  (2600  m2  g'1)  have  been  achieved  with  thermal 
activation.  After  polarization,  anions  and  cations  are  removed  from  the  electrolyte  by  the 
imposed  electric  field  and  electrosorbed  onto  the  carbon  aerogel.  The  solution  is  thus  separated 
into  two  streams,  concentrate  and  purified  water.  Based  upon  this  analysis,  it  is  concluded  that 
carbon  aerogel  CDI  may  be  an  energy-efficient  alternative  to  electrodialysis  and  reverse  osmosis 
for  the  desalination  of  brackish  water  (<  5000  ppm),  provided  that  cell  geometries  and  aerogel 
properties  are  carefully  tailored  for  such  applications.  The  intrinsic  energy  required  by  this 
process  is  approximately  QV/2,  where  Q  is  the  stored  electrical  charge  and  V  is  the  voltage 
between  the  electrodes,  plus  losses  due  to  parasitic  electrochemical  reactions,  electrical  resistance 
and  pressure  drop.  The  estimated  requirement  for  desalination  of  a  2000  ppm  feed  is  estimated 
to  be  -0.53-2.5  Wh  gal’1  (0.50-2.4  kJ  L’1),  depending  upon  voltage,  flow  rate,  cell  dimensions, 
carbon  aerogel  density,  recovery  ratio  and  other  parameters.  These  estimates  assume  that  50- 
70%  of  the  stored  electrical  energy  is  reclaimed  during  regeneration  (electrical  discharge).  The 
possibility  of  such  low  power  requirements  for  desalination  of  brackish  water  (BW),  as  well  as 
the  possibility  of  energy  storage  and  recovery,  may  make  this  process  attractive  for  such 
applications.  Though  the  intrinsic  energy  requirement  for  desalination  of  sea  water  (SW)  are  also 
relatively  low,  this  application  will  be  much  more  difficult.  Additional  work  will  be  required  to 
determine  the  suitability  of  carbon  aerogel  CDI  for  desalination  of  streams  that  contain  more  than 
5000  ppm  total  dissolved  solids  (TDS).  Applications  at  2000  ppm  will  require  the  construction 
of  electrochemical  cells  with  extremely  tight,  demanding  tolerances.  At  the  present  time,  the 
process  is  best  suited  for  streams  with  relatively  dilute  impurities,  as  recently  demonstrated 
during  a  field  test  at  LLNL  Treatment  Facility  C. 


DTIC  QUALITY  INSPECTED  I 
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Introduction 

Background.  Throughout  time  man  has  used  fresh  water  for  drinking,  industrial,  and  agricultural 
purposes,  and  has  settled  where  suitable  water  was  available,  as  stated  by  Summers  [1]. 
Progressive  industrialization  and  expansion  of  irrigation  agriculture  add  to  the  ever  increasing 
use  of  fresh  water.  The  rapid  increase  of  the  global  population  and  the  non-uniform  distribution 
of  fresh  water  around  the  world  has  motivated  research  into  the  development  of  various 
desalination  methods.  If  pure  water  could  be  economically  obtained  from  SW,  it  would  have  a 
dramatic  affect  on  our  future  standard  of  living,  as  competition  for  fresh  water  increases. 
President  John  F.  Kennedy  said,  “If  we  could  ever  competitively,  at  a  cheap  rate,  get  fresh  water 
from  salt  water,  ...  this  would  be  in  the  long-range  interests  of  humanity,  ...  (and)  would  dwarf 
any  other  scientific  accomplishment.”  Desalination  of  SW  and  BW  offers  great  potential  for 
increasing  the  availability  of  fresh  water.  The  concentration  of  salt  in  SW  is  typically  35,000 
ppm,  though  it  may  reach  50,000  ppm  in  some  areas.  The  concentration  of  salt  in  BW  rarely 
exceeds  10,000  ppm,  while  most  BWs  in  California  have  concentrations  between  800  and  3200 
ppm.  These  levels  must  be  reduced  to  less  than  500  ppm  for  drinking  water.  In  some  cases, 
lower  product  concentrations  are  required  for  agricultural  applications.  Until  now,  we  have  been 
paying  a  low  price  for  fresh  water,  however  the  cost  may  increase  dramatically  as  demand 
increases  and  we  are  forced  to  use  other  methods  to  avail  ourself  of  fresh  water. 

Desalination  Processes 

Desalination  Options.  Most  large-scale  desalination  processes  in  the  world  are  based  on 
variations  of  evaporation  and  distillation.  These  energy-intensive,  thermal  processes  require  heat 
that  is  derived  from  burning  fossil  or  nuclear  fuel.  Despite  the  huge  energy  requirements, 
thermal  processes  are  favored  for  SW  desalination  since  large  production  rates  are  possible  with 
reasonable  investments  in  capital  equipment.  Electrodialysis  (ED),  electrodialysis  reversal 
(EDR)  and  reverse  osmosis  (RO)  systems  are  more  energy  efficient,  but  require  expensive  and 
troublesome  membranes.  Membranes  are  plagued  by  chemical  degradation,  biological  fouling 
and  inorganic  scale  formation  problems,  and  must  be  periodically  replaced.  Membrane  processes 
have  been  used  for  desalination  of  SW,  but  are  more  often  considered  for  the  desalination  of 
reservoirs  of  brackish  water  (BW).  In  some  cases,  capacitive  deionization  (CDI)  with  carbon 
aerogel  electrodes  may  serve  as  an  energy-efficient  alternative  to  thermal  and  membrane 
desalination  processes.  Carbon  aerogel  is  an  ideal  electrode  material  because  of  its  low  electrical 
resistivity  (<  40  mQ  cm),  high  specific  surface  area  (2.0-5. 4xl06  ft2  lb'1  or  400-1100  m2  g’1), 
controllable  pore  size  distribution  (<  50  nm),  and  monolithic  structure.  Note  that  surface  areas  as 
high  as  1 ,3xl07  ft2  lb"1  (2600  m2  g"1)  have  been  achieved  with  activation. 

Thermal  Processes.  Thermal  processes  are  the  oldest  and  most  commonly  used  methods  of 
desalination  in  the  world.  The  thermal  efficiency  of  such  a  process  is  usually  called  the 
“performance  ratio”  or  “economy”  and  is  usually  defined  as  the  pounds  of  distillate  produced  per 
1000  Btu  of  heat  input.  Here,  the  total  required  energy  (electrical  and  thermal)  is  given  in  units 
of  Wh  gal 1  (kJ  L'1).  Evaporative  distillation  serves  as  the  major  source  of  fresh  water  in  the 
Middle  East.  Early  designs  are  of  the  submergedrtube,  multiple-effect  type.  In  such  systems, 
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steam  is  fed  through  tubes  submerged  in  a  stagnant  pool  of  brine,  thereby  causing  the  brine  to 
boil.  Evaporated  water  is  condensed  as  pure,  salt-free  water.  Salts  remain  in  the  brine  phase, 
which  becomes  more  concentrated.  Energy  recovery  is  accomplished  by  heat  transfer  from  the 
product  to  the  feed.  Scale  formation  on  heat  transfer  surfaces  lowers  the  thermal  efficiency  of 
the  process  and  is  exacerbated  by  high  operating  temperatures,  which  leads  to  localized 
supersaturation  and  precipitation.  Multistage  flash  evaporation  (MSF)  is  an  alternative  thermal 
process  that  was  developed  to  improve  process  efficiency  and  to  minimize  operating  problems 
associated  with  scale  formation.  Since  the  introduction  of  MSF  in  the  1960’s,  about  56%  of  the 
total  installed  desalination  plant  capacity  has  been  MSF  [2],  In  this  process,  salt  water  is 
introduced  into  a  low  pressure  chamber  where  it  undergoes  flash  evaporation.  Note  that  vapor 
can  be  produced  from  a  liquid  at  its  boiling  point  by  adding  heat  at  constant  pressure  and 
temperature  (boiling),  or  by  reducing  pressure  at  constant  temperature  (flashing).  The  amount  of 
heat  that  is  lost  from  the  process  is  minimized  by  using  several  counter-current  stages,  each  at  a 
lower  temperature  and  pressure  than  the  preceeding  stage.  Energy  recovery  is  improved  by 
adding  stages.  Since  heat  transfer  surfaces  in  MSF  processes  are  maintained  at  lower 
temperatures  than  those  in  submerged-tube  multiple-effect  systems,  problems  associated  with 
scale  formation  are  less.  A  typical  MSF  system  consumes  -320  Wh  gal'1  (304  kJ  L'1)  for  SW 
desalination,  compared  to  the  theoretical  minimum  given  here  of  -2.0  Wh  gal'1  (1.9  kJ  L'1). 
Note  that  a  simple  evaporator  without  energy  recovery  would  require  more  than  2720  Wh  gal'1 
(2587  kJ  L  ').  Energy  recovery  can  be  enhanced  by  using  mechanical  vapor  recompression 
(MVC).  Mechanical  compression  is  used  to  elevate  the  temperature  of  vapor  from  the  evaporator 
to  a  point  where  it  can  be  used  to  drive  the  evaporator.  The  hot  compressed  vapor  loses  its  latent 
heat  to  the  incoming  feed  and  condenses,  yielding  salt-free  product.  Since  the  temperature 
differential  between  the  hot  and  cold  side  of  the  heat  transfer  surface  is  only  7-9  °F  (4-5  °C),  the 
required  energy  and  rates  of  scale  formation  are  both  low.  A  typical  MVC  system  consumes 
-30-41  Wh  gal 1  (29-39  kJ  L  *).  See  Table  1  for  published  energy  requirements  for  SW 
desalination  processes  [2]. 

Reverse  Osmosis.  Reverse  osmosis  (RO)  is  a  separation  processes  which  uses  membranes  to 
selectively  remove  water  from  a  salt  solution.  In  conventional  osmosis,  water  or  solvent  flows 
through  a  semipermeable  membrane  from  a  less  concentrated  solution  to  an  area  of  higher 
concentration.  This  normal  osmotic  flow  can  be  reversed  by  applying  pressure  in  excess  of  the 
osmotic  pressure.  The  flow  through  the  semipermeable  membrane  is  governed  by  Equation  1 : 

=  A(Ap-  An)  [1] 

where  Fw  is  the  water  flux  (g  cm 2  s'1),  A  is  the  water  permeability  constant  (g  cm'2  s'1  atm'1),  Ap 
is  the  differential  pressure  applied  across  the  membrane  (atm),  and  An  is  the  osmotic  pressure 
differential  across  the  membrane  (atm)  [3a],  The  osmotic  pressure  is  a  function  of  salt 
concentration  and  is  represented  by  Equation  2: 

An  =  1.12(7’  +  273)^«ry  [2] 
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where  Arris  the  osmotic  pressure  (psi),  T  is  the  temperature  (°C),  and  mt  is  the  molality  of  the  i- 
th  ionic  or  non-ionic  constituent  [3b].  A  rule  of  thumb  that  works  well  with  natural  water,  based 
on  NaCl,  is  that  the  osmotic  pressure  increases  by  approximately  0.01  psi  (6.8x1  O'4  atm)  for  each 
ppm  (mg  L’1)  increase  in  salt  concentration.  Applying  this  rule  to  SW,  the  osmotic  pressure  is 
estimated  to  be  about  350  psi  (24  atm).  The  salt  flow  through  a  semipermeable  membrane  can  be 
expressed  by  Equation  3: 

Fx  =  B(Q-C2)  =  BAC  [3] 

where  B  is  the  salt  permeability  constant  (cm  s’1)  and  AC=CrC7  is  the  concentration  gradient 
across  the  membrane  (g  cm’ )  [3a],  RO  with  cellulose  acetate,  polyamide  or  other  polymer 
semipermeable  membranes  can  be  used  to  desalinate  BW  or  SW.  Usually,  membranes  are 
incorporated  into  plate-and-frame,  tubular,  or  spiral-wound  modular  units.  An  RO-based 
desalination  process  usually  consists  of  several  of  these  modules  connected  together  in  series- 
parallel  arrays.  The  processed  feedwater  is  separated  into  two  streams,  permeate  and  brine.  The 
permeate  passes  through  the  membranes  due  to  the  pressure  gradient  while  the  brine  is 
continuously  rejected.  Typical  systems  reject  brine  with  a  salt  concentration  of  two  to  four  times 
(2X  to  4X)  greater  than  that  of  the  feed.  Since  mechanical  energy  is  used  in  lieu  of  heat  energy, 
RO  systems  operate  at  ambient  temperature  and  avoid  some  of  the  problems  associated  with 
high-temperature  thermal  processes.  RO  systems  do  suffer  from  problems  such  as  concentration 
polarization,  chemical  degradation,  biological  fouling,  and  scaling  of  membranes.  Concentration 
polarization  at  the  brine-membrane  interface  increases  the  mass  transfer  resistance  and  requires 
additional  pressure  to  overcome.  Furthermore,  it  causes  localized  supersaturation  and 
precipitation,  which  exacerbates  scale  formation.  As  pores  in  the  membrane  become  blocked, 
the  flux  of  water  at  a  constant  applied  pressure  will  reduce.  In  order  to  avoid  such  problems, 
feed  water  must  be  pretreated.  RO  systems  are  generally  classified  as  sea  water  (800-1500  psi), 
standard  pressure  (400-650  psi),  low  pressure  (200-300  psi)  and  nanofiltration  (45-150  psi) 
[3a,3b],  Estimates  of  the  osmotic  pressure,  membrane  pressure,  and  minimum  possible  energy 
requirement  for  desalination  of  feed  streams  with  1000-35,000  ppm  TDS  are  found  in  Table  2. 
These  estimates  are  based  upon  heuristic  guidelines.  A  typical  RO  system  requires  -25-36  Wh 
gal  (24-34  kJ  L  )  for  desalination  of  SW,  depending  upon  the  use  of  energy  recovery,  whereas 
a  low-pressure  RO  system  requires  -8.5  Wh  gal'1  (8.1  kJ  L'1)  for  desalination  of  BW  with  1630 
ppm  TDS  [2,4].  Extremely  large  RO  plants  have  been  built  and  successfully  operated.  For 
example,  the  Yuma  Desalting  Plant  in  Yuma,  Arizona  is  the  largest  RO  facility  in  the  world  and 
is  operated  by  the  U.S.  Department  of  Interior,  Bureau  of  Reclamation.  Construction  of  this 
plant  was  authorized  by  the  Colorado  River  Basin  Salinity  Control  Act  (PL  93-320),  legislation 
passed  by  the  U.S.  Congress  in  1977.  The  objective  of  this  monumental  project  is  to  enable  the 
United  States  to  satisfy  conditions  of  a  treaty  with  Mexico  (Treaty  TS  994,  1944;  Ammended 
with  Minute  No.  242,  1977).  Under  that  agreement,  the  United  States  agreed  that  the  salinity  of 
water  delivered  to  Mexico  at  Morelos  Dam  would  not  exceed  the  salinity  level  at  the  Imperial 
Dam  by  more  than  1 15±39  ppm.  The  Imperial  Dam  is  27  miles  (43  km)  upstream  of  the  Morelos 
Dam.  The  Yuma  Desalting  Plant  uses  cellulose  acetate  RO  membranes  to  desalinate  as  much  as 
72.4xlOb  gal  d  (275xl06  L  d’1)  of  saline  drainage  water  from  farmlands  east  of  Yuma.  The  salt 
concentration  is  lowered  from  -3000  to  <  300  ppm  with  a  differential  pressure  of  -362  psi  (25 
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atm).  The  product  (-300  ppm)  and  unprocessed  drainage  water  (~3000  ppm)  are  blended  and 
sent  to  the  Colorado  River,  while  concentrated  brine  (-10,000  ppm)  is  rejected  and  sent  to  the 
Santa  Clara  Marsh  at  the  Gulf  of  California.  This  RO-based  desalination  plant  enables  the 
United  States  to  salvage  drainage  water  that  otherwise  would  be  too  saline  to  deliver  to  Mexico, 
thereby  saving  up  to  78,500  acre-feet  (97x1 06  mJ)  of  Colorado  River  water  per  year.  The  Yuma 
Desalting  Plant  is  not  being  operated  at  the  present  time  since  water  quality  requirements  can  be 
temporarily  met  without  it.  It  remains  in  a  stale  of  readiness. 

Electrodialysis.  In  a  conventional  electrodialysis  (ED)  cell,  several  pairs  of  anion  and  cation 
exchange  membranes  are  placed  between  a  pair  of  planar  electrodes,  one  serving  as  the  anode 
and  the  other  serving  as  the  cathode.  Isolated  compartments  are  formed  between  each  adjacent 
pair  of  anion  and  cation  exchange  membranes.  The  electric  field  imposed  by  polarization  of  the 
cell  causes  the  electromigration  of  ions,  with  negatively-charged  anions  moving  through  anion 
exchange  membranes  towards  the  positively-charged  anode,  and  positively-charged  cations 
moving  through  the  cation  exchange  membrane  towards  the  negatively-charged  cathode.  Ions 
are  forced  out  of  compartments  with  the  same  polarity  as  the  cell  (even-numbered)  and 
concentrated  in  compartments  having  opposite  polarity  (odd-numbered).  Relatively  pure  water  is 
withdrawn  from  even-numbered  compartments,  while  the  brine  phase  is  recycled  to  the  odd- 
numbered  compartments.  Electrolysis  causes  oxygen  evolution  at  the  anode  and  hydrogen 
evolution  at  the  cathode.  These  gaseous  reaction  products  must  be  continuously  removed  from 
the  cell.  Other  electrochemical  reactions  are  also  possible,  depending  upon  the  impurity  ions. 
Water  in  the  anode  compartment  becomes  acidic  due  to  the  accumulation  of  hydrogen  ions, 
while  water  in  the  cathode  compartment  becomes  alkaline  due  to  the  accumulation  of  hydroxyl 
ions.  Elevated  pH  in  the  cathode  compartment  promotes  precipitation  and  accelerates  scale 
formation,  and  makes  rinsing  of  the  electrodes  mandatory.  Extensive  pretreatment  of  the 
feedwater  is  also  required  to  reduce  membrane  degradation  and  fouling.  Process  costs  are 
directly  related  to  the  concentration  of  salt  in  the  feedwater,  ie..  the  amount  of  current  necessary 
to  produce  the  separation.  Electrodialysis  reversal  (EDR)  is  a  variant-  of  ED  with  periodic 
voltage  reversal  to  mitigate  scale  formation  and  fouling.  In  regard  to  Table  3,  it  was  assumed 
that  EDR  requires  about  2.0  Wh  gal"1  (1.9  kJ  L"1)  per  1000  ppm  TDS,  with  an  additional  2.5  Wh 
gal 1  (2.4  kJ  L'1)  for  pumping  and  an  allowance  of  5%  for  instrumentation  [3a].  Table  3  indicates 
that  EDR  requires  5.8  Wh  gal’1  (5.5  kJ  L'1)  to  desalinate  a  1500  ppm  BW  stream,  and  76.1  Wh 
gal  (72.4  kJ  L’1)  to  desalinate  SW.  Table  4  summarizes  the  estimated  energy  requirements  for 
two  types  of  BW  desalination  plants,  RO  and  EDR,  and  shows  that  EDR  requires  -7.7  Wh  gal"1 
(7.3  kJ  L"1)  for  desalination  of  a  feed  stream  with  1600  ppm  TDS  [4],  Other  references  indicate 
that  60-75  Wh  gal"1  (57-71  kJ  L"1)  is  required  by  EDR  for  desalination  of  SW  [5,6], 

Flow-Through  Capacitors  with  Activated  Carbon  Electrodes.  There  were  attempts  to  desalinate 
water  with  activated  carbon  electrodes  in  the  early  1960’s,  however,  limitations  of  the  electrode 
materials  available  at  that  time  prevented  further  development.  The  University  of  Oklahoma 
appears  to  have  been  the  first  group  to  use  flow-through  capacitors  with  activated  carbon 
electrodes  for  desalination,  with  seminal  publications  that  appeared  in  the  early  1960’s  [7,8]. 
Activated  carbon  powders  and  fibers  were  held  together  in  electrodes  by  a  variety  of  polymeric 
binders.  Johnson  et  al.  conducted  similar  studies  of  reversible  electrosorption  with  beds  of 
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activated  carbon  in  the  early  1970’s  [9-11].  Johnson’s  work  prompted  Newman  to  develop  a 
comprehensive  theoretical  model  for  the  capacitive  charging  of  porous  carbon  electrodes  [12]. 
Additional  work  was  done  in  Israel  several  years  later  and  published  in  the  1980’s  [13-15].  The 
work  in  Israel  viewed  the  process  as  a  chromatographic  separation.  Several  practical  problems 
were  encountered  with  these  conventional  activated  carbon  systems.  For  example,  the 
performance  (electrosorption  capacity)  of  activated  carbon  was  found  to  degrade  with  time.  In 
electrodes  made  with  polymeric  binders,  significant  fractions  of  the  activated  carbon  surface  was 
occluded.  Electrochemical  cells  that  used  flow-through  beds  of  activated  carbon  powder  as 
electrodes  required  membrane  separators  for  electrical  insulation  and  to  prevent  entrainment  of 
individual  particles  in  the  flow.  Furthermore,  flow  through  such  porous  media  is  characterized 
by  high  pressure  drop.  Process  efficiency  is  lowered  by  large  potential  drops  that  develop  in 
thick  electrodes,  which  include  packed  beds.  Finally,  early  activated-carbon  systems  had  no 
provisions  for  the  continuous  production  of  pure  water,  or  for  the  recovery  of  stored  electrical 
energy.  Eventually,  this  approach  to  desalination  was  abandoned. 

Capacitive  Deionization  with  Carbon  Aerogel  Electrodes 

Principle  of  Operation.  More  recently,  an  electrochemical  process  for  the  capacitive 
deionization  (CDI)  of  water  with  stacks  of  carbon  aerogel  composite  (CAC)  electrodes  has  been 
developed  by  Lawrence  Livermore  National  Laboratory  [16-24],  An  aqueous  solution  of  NaCl, 
Na2C03,  Na2S04  or  another  salt  is  passed  between  numerous  pairs  of  carbon  aerogel  electrodes, 
each  having  a  very  high  BET  surface  area  (2.0-5.4xl06  ft2  lb'1  or  400-1100  m2  g'1)  and 
exceptionally  low  electrically  resistivity  (<  40  mQ  cm).  Measured  electrosorption  capacities, 
extrapolated  to  0.03  N,  are  on  the  order  of  11.6xl0"3  lbNaC,  lbc'1  (20xl0'5  equiv  g'1).  After 
polarization,  ions  such  as  Na+  and  Cl"  are  removed  from  the  electrolyte  by  the  imposed  electric 
field  and  held  in  electric  double  layers  formed  at  the  surfaces  of  electrodes.  The  Gouy-Chapman 
Theory  or  the  Stem  Modification  of  that  theory  can  be  used  to  describe  the  effect  of  an  imposed 
electric  field  on  the  surface  charge  density  of  the  electric  double  layer  [25,26],  As  desired,  the 
effluent  from  the  cell  is  purified  water.  This  process  is  also  capable  of  simultaneously  removing 
a  variety  of  other  impurities.  For  example,  dissolved  heavy-metal  ions  can  be  removed  by 
reversible  electrosorption  [27]  or  by  electrodeposition  [28-31],  in  contrast  to  simple  double-layer 
charging.  In  the  specific  case  of  carbon  aerogel  CDI,  LLNL  has  investigated  a  broad  range  of 
cell  voltages  [16-24].  The  best  performance  (salt  removal)  was  achieved  at  -1.2-1. 3  V,  though 
operation  at  0.6  V  was  found  to  be  possible.  Successful  testing  was  done  at  solution 
conductivities  of  4,  59,  294,  588  and  5,882  ppm  (7,  100,  500,  1000  and  10,000  pS  cm"1).  In 
routine  single-pass  experiments,  more  than  99%  of  the  salt  was  removed  from  a  59  ppm  (100  pS 
cm  )  feed  stream.  After  the  carbon  aerogel  electrodes  became  saturated  with  salt,  breakthrough 
was  observed.  Electrodes  are  regenerated  by  electrical  discharge  prior  to  breakthrough  in  process 
applications,  which  allows  the  captured  salt  ions  to  be  released  into  a  relatively  small, 
concentrated  purge  stream. 

Carbon  Aerogel  Electrodes.  Carbon  aerogel  is  an  ideal  electrode  material  because  of  its  low 
electrical  resistivity  (<  40  mQ  cm),  high  specific  surface  area  (2.0-5.4xl06  ft2  lb"1  or  400-1100 
m  g"  ),  and  controllable  pore  size  distribution  (<  50  nm).  Note  that  thermal  and  chemical 
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activation  has  been  used  to  achieve  BET  surface  areas  of  1.3xl07  ft2  lb’1  (2600  m2  g" ' )  [32]. 
Resorcinol-formaldehyde  (RF)  aerogels  and  their  carbonized  derivatives  were  first  developed  by 
LLNL  [33-36],  Monolithic  sheets  of  this  material  were  made  by  infiltrating  a  resorcinol- 
formaldehyde  solution  into  a  porous  carbon  paper  (Textron  Specialty  Materials.  Lowell.  MA). 
curing  the  wetted  paper  between  glass  plates  in  a  closed  vessel,  and  then  pyrolizing  in  an  inen 
atmosphere.  The  exceptionally  high  conductivity  of  this  carbon  aerogel  composite  (CAC),  in 
contrast  to  loosely  bonded  carbon  powders  or  activated  carbon  fiber  cloths  (ACFCs),  is 
attributable  to  its  monolithic  structure  which  is  composed  of  interconnected,  covalently-bonded 
carbon  particles,  each  having  a  diameter  of  -12  nm.  In  contrast  to  electrodes  made  from 
activated  carbon  powders  and  fibers,  the  activation  energy  for  carrier  transport  in  carbon  aerogel 
is  relatively  small.  These  very  desirable  characteristics  have  enabled  LLNL  to  use  CAC  sheets 
as  electrodes  in  novel  supercapacitors  with  high  energy  density  and  high  power  density  [37], 
However,  these  energy-storage  devices  were  not  designed  to  permit  electrolyte  flow  and  required 
membranes  to  physically  separate  the  electrodes.  As  discussed  here,  LLNL  has  now  used  CAC 
electrodes  in  a  variety  of  configurations  to  remove  ionic  contaminants  from  water.  The  CAC 
used  to  generate  the  data  shown  in  this  publication  was  produced  at  LLNL.  as  well  as  a: 
GenCorp- Aerojet  Plant  in  Sacramento,  California.  During  this  first  commercial  production 
campaign,  it  was  demonstrated  that  a  large  quantity  of  high-quality  CAC  sheets  (4000  ft2  or  372 
m  )  could  be  mass  produced  at  an  apparent  cost  of  approximately  $50  per  square  foot.  It  may  be 
possible  to  reduce  the  cost  to  well  below  $  1-2  per  square  foot  in  the  future. 

Electrochemical  Cells.  Double-sided  electrodes  for  the  electrochemical  cell  were  made  by 
gluing  two  sheets  of  the  CAC  to  both  sides  of  a  titanium  plate  with  graphite-filled  epoxy.  The 
titanium  plate  served  as  both  a  current  collector  and  a  structural  support.  Each  sheet  of  CAC  was 
4  in  x  8  in  x  0.010  in  (10.16  cm  x  20.32  cm  x  0.0254  cm)  and  had  a  total  estimated  BET  surface 
area  of  ~3.4xl 04  ft2  (3.2xl07  cm2).  A  typical  stack  of  150  double-sided,  titanium-supported 
electrodes  had  a  total  estimated  BET  surface  area  of  approximately  — l.OxlO7  ft2  (9.4xl09  cmA. 
Electrolyte  flowed  through  the  stack  in  open  channels  formed  between  adjacent  electrodes.  An 
electrode  separation  of  -0.027  in  (0.069  cm)  was  maintained  in  these  early  cells  by  rubber 
compression  seals.  By  arranging  the  electrodes  so  that  orifices  alternate  from  one  side  of  the 
stack  to  the  other,  flow  from  the  bottom  of  the  stack  to  the  top  was  serpentine.  In  desalination 
applications,  the  titanium  will  have  to  be  replaced  with  a  less  expensive  alternative  and  smaller 
electrode  separations  (gaps)  will  be  required. 

Automated  Control  and  Potential-Swing  Operation.  A  prototypical  potential-swing  system  has 
been  developed  [24,38].  Ultimately,  such  synchronous  operation  will  be  essential  for  energv 
recovery,  as  well  as  continuous  production,  and  will  require  the  degree  of  automation  and 
sophistication  that  has  now  been  demonstrated.  This  prototype  can  produce  uninterrupted  flow's 
of  product  and  concentrate  by  operating  two  stacks  of  electrodes  in  parallel.  One  stack  purifies 
while  the  other  is  electrically  regenerated.  Flow  is  generated  by  a  programmable,  magneticallv- 
coupled,  screw  pump  with  a  304  stainless  steel  head.  All  lines  are  made  of  Teflon  and  had  a 
nominal  diameter  of  1/4  in  (0.635  cm).  The  cells  are  polarized  by  programmable  power  supplies 
that  have  a  voltage  range  of  0  to  12  V  or  a  current  range  of  0  to  60  A.  Sensors  were  placed  on 
the  inlet  and  outlet  lines  of  the  electrode  stack.  .  Electrical  conductivity  of  the  solution,  pH. 
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individual  ion  concentrations,  and  temperature  are  continuously  monitored.  The  system  is 
controlled  by  a  personal  computer.  A  single  AT-MIO-16DH  data  acquisition  board,  installed  in 
the  computer,  provides  the  interface  to  the  Input-Output  (I/O)  Signal  Subsystem.  The  I/O 
subsystem  consists  of  a  single  12-slot  Signal  Conditioning  Extension  Interface  (SCXI)  chassis. 
The  chassis  contains  seven  8-channel  analog-to-digital  (A/D)  modules  to  measure  flow,  level, 
pressure,  temperature,  pH,  and  conductivity;  two  16-channel  single-pole  double-throw  (SPDT) 
relay  modules  for  controlling  pumps  and  valves;  and  one  6-channel  digital-to-analog  (D/A) 
module  for  controlling  power  supply  voltages  and  pump  speed.  The  SCXI  hardware  provides 
multiplexing,  filtering,  isolation,  and  amplification  for  the  process  signals.  The  operating  system 
is  DOS  v6.22  running  Windows  v3.1.  LabVIEW  v3.1  software,  running  under  Windows,  is  used 
for  data  acquisition  and  control.  LabVIEW  is  a  graphical  programming  environment  which 
provides  integrated  tools  for  acquisition,  control,  analysis,  and  presentation,  as  well  as 
connectivity  to  serial,  parallel,  voltage,  current  loop,  RTD,  thermistor,  and  relay  communication 
interfaces.  Data  acquisition  and  control  software  is  optimized  and  converted  to  compiled  run¬ 
time  code.  Operator  input  is  via  a  mouse  and  keyboard. 

Surface  Area  and  Electrosorption  Capacity 

Interpretation  of  BET  Surface  Areas.  It  is  noteworthy  that  activated  carbon  powders  with 
Brunauer-Emmet-Teller  (BET)  surface  areas  as  high  as  1 .5x10'  ft2  lb"1  (3000  m2  g"1)  are  readily 
available.  However,  much  of  the  surface  area  in  such  materials  is  located  inside  pores  having 
diameters  less  than  1  nm.  It  is  believed  that  the  electrochemically  active  area  is  only  a  fraction  of 
the  BET  surface  area.  BET  analyses  are  probably  misleading  since  gas  molecules  can  penetrate 
much  smaller  pores  than  a  typical  electrolyte.  For  example,  the  bond  length  of  N2  is  only  0.1 
nm.  It  is  very  doubtful  that  this  level  of  porosity  contributes  to  electrochemical  double  layer 
formation  since  electrolyte  penetration  and  double  layer  formation  are  questionable  on  this  scale. 
From  the  Gouy-Chapman  theory,  as  well  as  the  Stern  modification  of  that  theory,  it  is  believed 
that  a  fully-developed  electric  double  layer  on  a  planar  electrode  with  no  detrimental  shielding 
effects  would  require  much  greater  distances  for  full  development.  In  the  case  of  a  1:1 
electrolyte  in  water  at  25°C,  the  characteristic  thickness  of  the  diffuse  layer  ranges  from  1  nm  at 
a  concentration  of  0.1  M  to  30  nm  at  10"4  M  [26,27].  Though  BET  surface  areas  are  frequently 
quoted,  the  electrosorption  capacity,  given  as  equivalents  adsorbed  per  gram  of  carbon  is 
believed  to  be  a  more  relevant  measure  of  electrode  performance  [7,8]. 

From  the  Gouy-Chapman  theory  developed  for  simple  planar  electrodes,  one  might  expect  the 
surface  charge  density  to  have  a  square  root  dependence  on  electrolyte  concentration  [25,26].  In 
the  case  of  dilute  aqueous  solutions  at  25°C,  the  following  expression  should  be  obeyed  by  both 
anodes  and  cathodes: 

aGC  =  1 1.7 Vc  sinh(19.5  z<j>0)  [4] 

where  aQC  is  the  surface  charge  density  (p.C  cm"  ),  C  is  the  electrolyte  concentration  (mol  L"  ),  z 
is  the  ionic  charge,  and  <p0  is  the  electrode  potential  (mV).  However,  double  layer  formation  on 
carbon  electrodes  is  much  more  complicated.  For  example,  cathode  capacities  are  much  higher 
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than  corresponding  anode  capacities  [7,8].  This  is  attributed  to  the  cation  affinity  of  carbonyl 
groups  on  the  electrode  surface.  Therefore,  the  electrosorption  capacity  of  an  electrochemical 
cell  built  with  porous  carbon  electrodes  is  probably  limited  by  its  ability  to  accommodate  anions. 

Recent  data  from  experiments  with  NaCl  are  summarized  in  Tables  5a  and  5b,  and  illustrated  by 
Figures  1  and  2.  Figure  1  shows  the  electrosorption  of  Na+  and  Cl  on  54  double-sided  CAC 
electrodes,  each  having  a  width  of  4  in  and  a  length  of  8  in  (10.16  cm  x  20.32  cm).  The  quantity 
of  electrosorbed  NaCl  (Y)  was  calculated  front  the  minimum  in  electrolyte  concentration  (X]), 
reached  after  polarizing  the  electrodes  for  about  55  minutes.  The  stack  of  electrodes  was 
maintained  at  a  cell  voltage  of  1.2  V  (X2)  while  0.93  gal  (3.5  L)  of  electrolyte  was  continuously 
recycled  at  a  rate  of  1.59  gal  h'1  (X3,  100  ml  min'1).  Results  for  59,  278  and  588  ppm  (100,  500 
and  1000  pS  cm'1)  NaCl  solutions  are  compared.  Figure  2  shows  electrosorption  on  100  double¬ 
sided  CAC  electrodes  with  the  same  dimensions  as  those  represented  by  Figure  1 . 
Approximately  0.79  gal  (3  L)  of  59  ppm  (100  pS  cm'1)  NaCl  solution  was  continuously  recycled 
at  a  rate  of  1.59  gal  h'1  (100  ml  min'1).  Results  for  0.6  and  1.2  V  are  compared.  Multiple 
variable  regression  analysis  of  this  data  was  used  to  determine  the  dependence  of  the  anion 
capacity,  Y,  on  equilibium  electrolyte  concentration,  X,,  and  cell  voltage,  X2.  The  effects  of 
flow  rate  were  not  taken  into  consideration.  The  results  are  summarized  by  Equation  5: 


Y  =  e'3959X?2'mX;-5m  [5] 

where  Y  has  the  units  of  1 0'5  equivalents  per  gram  of  carbon  aerogel  (based  on  the  anode  mass 
alone),  X,  has  the  units  of  10"3  equivalents  per  liter  of  solution,  and  X2  is  given  in  volts.  The 
agreement  between  the  correlation  and  experimental  data  is  reasonably  good,  as  illustrated  by 
Figure  3.  This  is  reflected  in  the  calculated  multiple-variable  regression  coefficient,  R  ,  which  is 
0.966.  A  second  regression  analysis  was  done,  which  also  accounted  for  the  effects  of  flow  rate, 
X3.  The  results  are  summarized  by  Equation  6: 


Y  = 


7.6730  v-0-27484  v2.7024 

e  A,  X2 


X- 


-0.27814 


[6] 


where  X3  is  given  in  ml  min'1.  Here  too  the  agreement  between  the  correlation  and  experimental 
data  is  good,  as  illustrated  by  Figure  4.  The  calculated  multiple- variable  regression  coefficient, 
R2,  is  0.942.  Since  Y  and  X,  are  directly  proportional  to  a  and  C,  respectively,  one  could 
reasonably  expect: 


Y  cc  AT” 5  [7] 

It  is  believed  that  the  failure  to  observe  square-root  concentration  dependence  is  due  to  self¬ 
shielding  effects  experienced  by  the  porous  carbon  electrodes. 

Breakthrough  in  a  Stack  of  Carbon  Aerogel  Electrodes 

The  total  area  of  carbon  aerogel  sheet  required  to  desalinate  a  given  feed  stream  (S)  can  be 
estimated  with  the  following  equation: 
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S«0(C,-C,)-S-  [8] 

<7 

where  0  is  the  volumetric  flow  rate,  CF  is  the  concentration  of  salt  in  the  feed  stream,  CP  is  the 
concentration  of  salt  in  the  product  stream,  a  is  the  area-specific  electrosorption  capacity  of  the 
carbon  aerogel  at  CF,  and  rB  is  the  breakthrough  time.  Equation  8  assumes  that  the  concentration 
shock  wave  is  an  abrupt  step  function.  Note  that  the  area-specific  electrosorption  capacity  is: 

[9] 

v/here  FWsall  is  the  formula  weight  of  the  salt,  paerogei  is  the  density  of  the  carbon  aerogel  sheet, 
<5 'aer0gei  is  the  thickness  of  the  carbon  aerogel  sheet,  and  Y  is  the  electrosorption  capacity, 
calculated  with  an  empirical  expression  established  by  linear  regression  analysis  of  actual  data. 
As  shown  in  Tables  6a  and  7a,  the  density  of  CAC  sheets  was  assumed  to  be  0.413  g  cm°  in 
Cases  1  through  5,  and  0.800  g  cm"3  in  Case  6.  Then,  for  example,  in  Cases  1  through  5: 

P aerogel  S  aerogel  ~  0.0  1  05  g  Cnf1  [10] 

It  should  be  noted  that  significant  variability  in  the  density  of  CAC  sheets  from  production 
campaigns  has  been  observed.  Another  important  time  constant  is  the  liquid-phase  residence 
time,  rR,  can  be  calculated  from  the  channel  volume,  V,  and  the  volumetric  flow  rate,  0. 


if  plug  flow  is  assumed.  The  liquid  volume  in  the  channel  is  assumed  to  be: 

r«*L(h  +  WS^)  [12] 

where  w  is  the  channel  (electrode)  width,  L  is  the  channel  length,  h  is  the  electrode  separation 
and  <9  is  the  void  fraction  in  the  carbon  aerogel  sheet.  Note  the  relationship  between  S,  L,  and  w. 


The  void  fraction  is  simply: 


q  _  j  Paemgel 
P  carbon 
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where  pcarbon  is  assumed  to  be  2.7  g  cm'3  (values  as  low  as  2.25  g  cm  3  can  be  found  in 
handbooks).  Obviously,  for  the  system  to  work  properly,  the  concentration  shock  wave  must 
move  through  the  stack  more  slowly  than  the  liquid: 


*7,  ^  tr 


[15] 


The  recovery  ratio,  R ,  is  defined  as  the  fraction  of  the  feed  recovered  as  desalinated  product  and 
can  be  expressed  in  terms  of  rB  and  tR. 


R  =  [16] 

tb  +  tr 

This  assumes  that  the  flow  rate  during  deionization  and  regeneration  are  equivalent.  In  an  ideal 
system,  the  maximum  possible  recovery  ratio  will  be  about  50%  at  the  point  where  the  two  time 
constants  are  identical.  As  rB  becomes  greater  than  zR,  the  recovery  ratio  will  increase.  This 
criteria  can  be  used  to  select  a  consistent  electrode  separation,  which  will  dictate  the  pressure 
drop  in  the  system. 


h  < 


1  -R  a 
R  (C,-C,,) 


-2  S 


aerogel 


(1  P aerogel  /  P  carbon  ) 


[17] 


As  the  concentration  in  the  feed  stream  increases,  the  anode-cathode  separation  must  become 
smaller.  It  is  important  to  note  combinations  of  a  and  CrCP  exist  that  make  it  impossible  to  find 
a  reasonable  separation  ( h  would  be  zero  or  negative).  Alternatively,  the  minimum  possible 
product  concentration  can  be  determined  for  a  fixed  electrode  separation. 


CP>CF- 


1  -R _ a _ 

R  k  +  2  3 acn)gei(  1  P aerogel  P carbon  ) 


[18] 


Clearly,  the  value  of  Cp  will  have  to  be  greater  than  zero,  a  requirement  that  places  constraints  on 

C/r?  R~ 5  o*,  hj  & 'aerogel  P aerogel'  Equations  17  and  18,  the  factor  2  Saer0gsi(l  ~Paeroge / Pcarbon) 

assumes  that  the  CAC  sheets  used  to  construct  anodes  and  cathodes  have  identical  theknesses. 
Since  the  process  may  be  anion  limited  [8],  one  might  be  able  to  use  much  thinner  CAC  sheets 
for  the  cathodes  than  for  the  anodes,  perhaps  with  Scalhode  ~  0.2Sanode.  In  such  a  case,  the  factor 
becomes  l-28aerogel(l-paeroge/pcarboJ,  where  Saeroge,  is  the  thickness  of  the  CAC  sheets  used  to 
make  the  anodes.  By  using  the  thinnest  possible  cathode,  larger  electrode  gaps  can  be  tolerated 
and  lower  product  concentrations  can  be  achieved. 


Minimum  Energy  Required  for  Separation 

The  minimum  work  of  separation  of  a  feed  mixture  into  impure  products  at  constant  temperature 
and  pressure  can  be  computed  with  Equation  19  [39]: 
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Kin  ,r  =~Rt{  Z  V  ZX/V  ln(r  ,7X/7  ) 

V  j  <  j  2 


where  $  is  the  molar  fraction  of  feed  entering  product  xy,  is  the  mole  fraction  of  component  j  in 
product  i,  and  is  the  activity  coefficient  of  component  j  in  product  /.  For  a  binary  mixture, 
Equation  19  becomes: 


W .  T  = 

min./ 


XA\  XA2 


(XAF  ~xaj)  XA\  lny  ^::-— +  (1  -  x^,)ln ^ 

L  Y  AiXAi  \-XAX 


KlXA!  +(1-^  )lnLidL 


+  (x^j  XA2 


The  Debye-Huckel  model  [40]  can  be  used  to  estimate  the  activity  coefficient  of  component  j: 


z]a'4f 


In  Yj  = - - - t= 

J  1  +  B'ay[T 


where  a  is  the  exclusion  distance  of  a  central  ion  located  at  the  origin,  z-  is  the  valence  of  species 
j,  I'  is  the  molar  ionic  strength,  and  a’  and  B  ’  are  the  Debye-Huckel  parameters.  The  molar  ionic 
strength  is: 


where  cy  is  the  molar  concentration  of  species  j.  The  Debye  Huckel  parameters  are  defined  as: 


,  Fe  n 

a  =  — - B 

8  7ieRT 


where  R  is  the  universal  gas  constant,  T  is  the  absolute  temperature,  F  is  Faraday’s  constant,  e  is 
the  electronic  charge,  and  £  is  the  permittivity.  As  the  ionic  strength  goes  to  zero,  the  Debye- 
Huckel  limiting  law  becomes: 


In  yj*-z)a.'JT 
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Equations  20  through  25  were  used  to  calculate  the  minimum  energy  requirements  for 
desalination  as  a  function  of  NaCl  concentration.  Results  are  summarized  in  Tables  6e,  6f  and 
7f.  The  minimum  energy  required  for  separation  of  a  2000  ppm  NaCl  stream  into  a  100  ppm 
product  and  a  concentrate  with  R=0.5  is  -0.095  Wh  gal"1  (0.090  kJ  L  '),  whereas  the  minimum 
energy  required  for  separation  of  a  35,000  ppm  NaCl  solution  (surrogate  SW)  into  a  100  ppm 
product  and  a  concentrate  with  R=0.5  is  -2.0  Wh  gal 1  (1.9  kJ  L  '). 

Electrical  Energy  Requirement 

The  energy  stored  in  the  electrochemical  cell  during  charging  (deionization)  can  be  estimated 
from. the  expression  for  an  ideal  capacitor: 

i I^=\QV  126] 

where  Q  is  the  stored  electrical  charge  and  V  is  the  cell  voltage.  Stored  energy  is  available  for 
recovery  during  discharge  (regeneration).  Equation  16  accounts  for  the  energy  consumed  by 
parasitic  electrochemical  reactions: 


U^  =  v\l„(t)dl  P7] 

0 

where  Ip(t)  is  the  instantaneous  current  due  to  parasitic  charge-transfer  reactions  and  l  is  the  time 
required  for  deionization.  Ohmic  losses  must  also  be  accounted  for: 


/ 

R(0[0(0  +  /„«)f<*  [28] 

0 

where  Ic(t)  is  the  instantaneous  current  due  to  deionization  and  R(t)  is  the  instantaneous  electrical 
resistance.  The  ohmic  loss  associated  with  the  parasitic  current  is  implicit  in  Equation  27,  and  is 
included  explicitly  in  Equation  28.  Consequently,  when  these  two  contributions  are  summed  to 
calculate  the  total  energy  consumption,  the  estimate  will  be  slightly  higher  than  the  actual  value 
(conservative).  Note  that  Ip(l),  Ic(t)  and  R(t)  are  spatially-averaged  quantities,  representing  the 
entire  stack  of  electrodes.  At  this  point,  a  reasonable  approximation  for  an  order-of-magnitude 
esimate  can  be  derived  by  treating  the  currents  and  the  resistance  as  a  time-average  quantities. 


U 


parasitic 


V  Ip  At 


[29] 


[30] 
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The  energies  per  unit  volume  of  product,  Ws,ored,  Wparasilic,  and  Wohmic,  are  calculated  by  dividing 
Uslored,  Uparasilic  and  Uohmic,  respectively,  by  the  volume  of  water  processed  per  unit  time: 


W„ 


Mo  red 


^ stored 

©At 


[31] 


W. 


parasitic 


U 


parasitic 


©At 


[32] 


l U  —  ^  ohmic 

r  ohmic 


©At 


[33] 


where  0  is  the  volumetric  flow  rate  of  product  during  the  interv  al  At. 

Parasitic  Current:  Limited  by  Diffusion  and  Fixed  Reactant  Mass 

The  observed  parasitic  current  may  be  due  to  a  fixed  mass  ot  electroactive  reactant  introduced 
with  the  feed.  One  possibility  might  be  the  reduction  of  dissolved  oxygen: 

02  +  2 H20  +  4e~  dcc'wJc  ■>  40H~  [34] 

which  has  a  standard  reduction  potential  of  0.401  V  measured  against  a  standard  hydrogen 
electrode  (SHE).  The  depletion  of  such  a  reactant  is  assumed  to  be  limited  by  mass  transport  to 
the  electrode-electrolyte  interface.  For  simplicity,  we  assume  that  this  interface  is  the  outer 
surface  of  the  monolithic  sheet  of  carbon  aerogel.  The  diffusive  flux  of  such  non-ionic  species  to 
the  surface  of  the  electrode  is: 


J  =  -D 


dC 

dy 


kAC 


[35] 


where  D  is  the  diffusivity  (~10’5  cm2  s"1),  C  is  reactant  concentration  in  the  bulk  electrolyte  (mol 
cm'3),  k  is  the  mass  transfer  coefficient  (cm  s"1),  and  y  is  the  direction  normal  to  the  carbon 
aerogel  electrode.  The  maximum  possible  flux  coincides  with  a  zero  surface  concentration  (AC 
=  Q.  Note  that  the  mass  transfer  coefficient  can  be  calculated  from  a  dimensionless  correlation 
such  as  the  Sherwood  equation: 


k°e  =  0.023 

'  DcuhP' 

0.8 

m) 

D 

V  gcp  9 

l  pD) 

The  term  on  the  LHS  of  tire  equation  is  known  as  the  Sherwood  number  (Sh),  the  first  term  on 
the  RHS  is  known  as  the  Reynolds  number  (Re),  and  the  second  term  on  the  RHS  is  the  Schmidt 
number  (Sc).  A  mass  balance  on  a  differential  length  of  the  electrochemical  cell  channel  yields: 
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dC 

dx 


[37] 


where  ©  is  the  volumetric  flow  rate.  Integration  of  this  expression  yields  the  reactant 
concentration  as  a  function  of  channel  length. 


C  =  CF  exp 


wkx 


[38] 


where  CF  is  the  concentration  of  the  reactant  in  the  feed,  w  is  the  electrode  (channel)  width,  and  x 
is  the  channel  length.  The  current  density  is  then  defined  in  terms  of  this  diffusive  flux: 

i-nFJ~nFkAC  [39] 


where  n  is  the  number  of  electrons  involved  in  the  electrochemical  reaction  and  F  is  Faradays 
constant.  The  limiting  current  is  calculated  from  the  maximum  possible  flux,  the  point  where  the 
surface  concentration  is  zero  (AC  =  C). 

h.  ~  nFJmsx  ~  nFkC  [40] 

Note  that  this  quantity  is  a  function  of  position  along  the  channel. 


iL(x)  =  nFkCF  exp 


wkx 

~e~ 


[41] 


The  total  parasitic  current  due  to  this  diffusion-limited  electrochemical  reaction  is  calculated  with 
the  following  integral: 


Ip  =  w\h(x)dx 

o 

The  integration  yields: 


[42] 


f 

wkx 

1-exp 

— 

V 

©  J \) 

[43] 


This  can  be  used  as  the  basis  of  scaling  a  measured  parasitic  current  to  any  channel  length: 
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1  -  exp 

1 

wkx. 

1 

© 

1  -  exp 

wkx2 

0 

The  parasitic  current  can  be  expressed  in  terms  of  the  total  number  of  stacks,  N. 


W'.A' 


I  p,n= i 


1-exp 


1  -  exp 


wk  Ax  N 
O 

wk  Ax 
0 


[44] 


[45] 


where  Ax  is  the  channel  length  of  a  single  stack.  The  maximum  possible  value  of  Ip  is  found  by 
taking  the  limit  as  N  approaches  infinity. 


lim  I r  A,  —  Ir  N= ! 


1 


A'— 


1  -  exp| 


vi'  k  Ax 
0 


[46] 


An  order-of-magnitude  estimate  of  the  parasitic  current  as  a  function  of  total  channel  length  or 
the  number  of  stacks  can  be  made.  For  example,  assume  that  h  >  0.069  cm  >  d,  w  -10.16  cm.  Ax 
-3000  cm,  D  -10"5  cm2  s’1  and  0  -100  cm3  min'1.  In  this  case,  the  parasitic  current  for  an 
infinite  number  of  stacks  is  only  5%  greater  than  the  parasitic  current  for  a  single  stack.  In  this 
case,  there  is  no  energy  penalty  for  increasing  the  channel  length. 

Parasitic  Current:  In  Situ  Generation  of  Reactants  Governed  by  Cell  Voltage 

The  observed  parasitic  current  might  also  be  due  to  very  low  levels  of  electrolysis  at  both 
electrodes.  In  this  scenario,  oxygen  generated  at  the  anode  would  dissolve  in  the  water.  After 
being  tranported  to  the  cathode,  it  would  undergo  reduction.  Hydroxyl  ions  produced  at  the 
cathode  during  oxygen  reduction  would  then  be  transported  to  the  anode,  where  oxidation  would 
provide  additional  oxygen.  Thus,  a  continuous  cycle  would  be  established  and  no  decline  in  the 
parasitic  current  density  over  the  length  of  the  channel  would  be  observed.  In  such  a  case,  it 
could  be  assumed  that  the  parasitic  current  density,  i  (pA  cm'2),  would  obey  the  following 
empirical  expression: 

ln(z)  =  A  +  BAV  [47] 

where  A  and  B  are  temperature-dependent  empirical  constants,  and  AV  (V)  is  the  cell  voltage. 
Some  preliminary  data  suggests  that: 

ln(r)=  1.221  +  1.5069  A  V  [48] 
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Based  on  this  empirical  equation,  the  parasitic  current  density  would  be  expected  to  be  about 
18.7  pA  cm'2  at  1.2  V,  and  7.6  pA  cm'2  at  0.6  V. 

Mechanical  Energy  Requirement 

The  mechanical  energy  balance  for  an  incompressible  fluid  is  given  as  Equation  49  [41-43]: 


Au[ 
2  gc 


£  A  p  , 

+  — A  z  + - +  Iw , 

gc  P 


+  77  W 

I  pump  pump 


=  0 


[49] 


where  ub  is  the  change  in  average  velocity  of  the  electrolyte  in  the  open  channel  between  pairs  of 
CAC  electrodes,  g  is  the  gravitational  acceleration,  gc  is  the  gravitational  constant,  z  is  the 
elevation,  p  is  the  pressure,  p  is  the  electrolyte  density,  Iwj  is  the  lost  work  due  to  friction,  rjpump 
is  the  efficiency  of  the  pump,  and  Wpump  is  the  energy  required  by  the  pump.  The  lost  work  due 
to  friction,  based  upon  the  Fanning  friction  factor,  is: 


2  fLuj 

gcDc 


[50] 


where  De  is  the  equivalent  diameter,  defined  as  four  times  (4X)  the  hydraulic  radius,  rh.  The 
hydraulic  radius  is  defined  as  the  cross-sectional  area  of  the  channel  divided  by  the  wetted 
perimeter: 


wx  h 
2  (w  +  h) 


[51] 


The  electrode  width  and  separation  are  w  and  h,  respectively.  The  Fanning  friction  factor  is  then 
defined  as: 


/ 


*  gc  Dc  (-A p) 

2  u]  p  L 


[52] 


Alternatively,  the  lost  work  due  to  friction,  based  upon  the  Moody  friction  factor,  is: 


‘  2  g,D, 


[53] 


The  Moody  friction  factor  is  defined  as: 


/ 


2gcDe(-Ap) 

u]pL 


[54] 
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Early  work  done  with  pipes  indicates  that  [41]: 


The  Reynolds  number  is  defined  as: 


Re  = 


E>eUbP 

ScP 


[55] 


[56] 


where  ju  is  the  viscosity.  A  regression  analysis  of  experimental  pressure-drop  data  indicates  that 
the  Moody  friction  factor  for  a  stack  of  150  double-sided  carbon  aerogel  electrodes  is: 


f  = 


14.3766 

Re0'72917 


[57] 


The  y  variance  was  0.206,  the  x-y  variance  was  0.01 1,  and  the  regression  coefficient  was  0.953, 
all  indicating  a  good  fit  to  the  data.  The  regression  equation  is  compared  to  the  actual  data  in 
Figure  5.  The  difference  between  Equations  55  and  57  is  probably  due  to  slippage  of  fluid  at  the 
boundary  between  the  aerogel  and  the  open  channel  where  the  water  is  assumed  to  flow.  More 
slippage  would  be  expected  at  a  fluid-filled  porous  wall  than  at  a  solid  wall,  and  would  lower  the 
effective  friction  factor.  Note  that  the  sheets  of  carbon  aerogel  were  4  in  x  8  in  x  0.010  in  (10.2 
cm  x  20.3  cm  x  0.025  cm);  the  epoxy  film  used  to  glue  the  carbon  aerogel  to  the  titanium 
substrate  was  approximately  0.0025  in  (0.0064  cm)  thick;  and  the  gaskets  used  to  separate  the 
titanium  plates  were  0.062  in  (0.16  cm)  thick  without  compression  and  about  0.052  in  (0.13  cm) 
thick  with  compression.  The  gap  separating  adjacent  sheets  of  carbon  aerogel  were 
approximately  0.027  in  (0.069  cm).  The  corresponding  hydraulic  radius  was  0.0134  in  (0.0341 
cm)  and  the  equivalent  diameter  was  0.0536  in  (0. 1 36  cm).  The  energy  requirement  for  pumping 
is  approximately  7.2xl0‘3  Wh  gaf 1  psf1. 

Total  Energy  Requirement 

The  total  energy  requirement  per  unit  volume  of  product  is  then  calculated  as: 

W'o,a>  =0-*)  Kand  +  ^ parasitic  +  Kittttic  +  ~Wp^P  [58] 

K 

where  e  is  the  fraction  of  the  stored  electrical  energy  that  is  recovered  during  discharge,  or 
regeneration  of  the  stack.  The  other  terms  are  defined  by  Equations  16,  31,  32,  33  and  49, 
respectively.  An  energy  recovery  of  approximately  50%  {s  ~  0.5)  should  be  relatively  easy  to 
achieve  since  no  voltage  conversion  should  be  required.  With  voltage  conversion,  approximately 
70%  ( s  ~  0.7)  should  be  possible.  As  already  discussed  in  regard  to  Equations  27  and  28,  the 
ohmic  loss  associated  with  the  parasitic  current  are  accounted  for  twice,  which  leads  to  a 
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conservative  estimate.  However,  the  result  is  inconsequental  since  the  ohmic  losses  are 
relatively  insignificant. 

Results  of  Energy  Analysis 

Theoretical  Minimum.  As  previously  discussed,  the  minimum  theoretical  work  required  by  an 
isothermal  process  to  separate  SW  or  BW  into  product  and  concentrate  streams  can  be  calculated 
with  equations  19  through  25.  If  SW  and  BW  are  assumed  to  be  NaCl  solutions,  Equation  20 
can  be  used  (as  was  done  here).  This  calculation  requires  that  the  activity  coefficient  of  NaCl  be 
estimated  with  the  Debye-Huckel  Theory.  Assuming  a  recovery  ratio  of  0.5,  the  theoretical 
minimum  work  for  separation  of  SW  (35,000  ppm  NaCl)  or  BW  (2000  ppm  NaCl)  into  product 
(100  ppm  NaCl)  and  concentrate  streams  are  ~2.0  Wh  gal'1  (1.9  kJ  L'1)  or  -0.10  Wh  gal'1  (0.09 
kJ  L’1),  respectively.  All  estimates  of  the  required  energy  for  carbon  aerogel  CDI  under  these 
conditions  should  be  greater  than  these  limiting  values. 

Estimated  Energy  Requirements  for  Carbon  Aerogel  CDI.  Cases  1  through  6  in  Tables  6a 
through  6f,  which  are  illustrated  by  Figures  6a  through  6f.  are  scenarios  where  the  minimum 
possible  product  concentration  was  fixed  at  100  ppm,  the  recovery  ratio  was  fixed  at  0.5  (except 
Case  6),  and  the  maximum  possible  electrode  gap  was  calculated.  Case  6  assumed  a  recovery 
ratio  of  0.85.  Cases  1  through  6  in  Tables  7a  through  7i,  illustrated  by  Figures  7a  through  7f,  are 
scenarios  where  the  electrode  gap  was  fixed,  which  enabled  calculation  of  the  maximum  possible 
recovery  ratio,  R,  and  the  minimum  possible  product  concentration,  Cp.  The  “Energy 
Requirement  -  No  Depletion”  assumes  that  the  electroactive  species  responsible  for  the  parasitic 
current  is  generated  in  situ  and  not  depleted,  whereas  the  “Energy  Requirement  -  With 
Depletion”  assumes  that  the  electroactive  species  is  introduced  in  the  feed  (dissolved  oxygen, 
etc.)  and  is  depleted  by  the  parasitic  electrode  reaction.  By  treating  the  channel  like  a  plug-flow 
reactor,  depletion  factors  can  be  calculated  (Equations  37  through  46],  Note  that  some 
reasonably  attractive  cases  for  brackish  water  were  identified,  but  no  good  combination  of 
parameters  for  sea  water  was  found.  In  particular.  Cases  3  and  6  in  Tables  7a  through  7i  appear 
to  have  relatively  low  energy  requirements,  as  well  as  reasonable  recovery  ratios.  As  specified  in 
Case  6,  it  will  be  necessary  to  use  thicker,  more  dense  CAC  sheets  {Saerogel  -  0.127  cm,  paeroge,~ 
0.8  g  cm"3)  than  currently  available  in  prototypes  to  achieve  high  recovery  ratios  ( R  -  0.863- 
0.890).  Such  a  system,  with  an  electrode  separation  of  about  13.8x1 0'3  in  (0.035  cm),  should  be 
able  to  desalinate  BW  (1500-2000  ppm)  with  an  energy  requirement  of  ~0.53- 1.06  Wh  gal"1 
(0.50-1.01  kJ  L"1).  This  case  is  compared  to  RO  and  EDR  in  Table  8. 

Comparison  of  Competing  Technologies  for  Desalination.  Competiting  technologies  for  SW 
desalination  include  reverse  osmosis  (RO),  mechanical  vapor  recompression  (MVC), 
electrodialysis  (ED),  multistage  flash  (MSF),  multiple  effect  distillation  with  mechanical  vapor 
recompression  (MED-MVC)  and  multiple  effect  distillation  with  thermal  vapor  recompression 
(MED-TVC).  Published  data  indicates  that  RO  consumes  -25-36  Wh  gal"1  (24-34  kJ  L"1),  MVC 
consumes  -30-41  Wh  gal'1  (29-39  kJ  L"1),  ED  consumes  -60-75  Wh  gal'1  (57-71  kJ  L"1),  MED- 
TVC  consumes  -215-315  Wh  gal'1  (205-300  kJ  L'1),  MSF  consumes  -320  Wh  gal'1  (304  kJ  L"1), 
and  simple  evaporation  without  energy  recovery  consumes  -2720  Wh  gal'1  (2586  kJ  L"1).  In  the 
case  of  RO,  lower  energy  requirements  can  be  achieved  through  the  use  of  energy  recovery 
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turbines.  The  intrinsic  energy  requirement  for  SW  desalination  with  carbon  aerogel  CDI  is  QV/2 , 
which  is  -36-37  Wh  gal'1  (34-35  kJ  L'1).  An  energy  requirement  less  than  that  needed  for  RO 
can  be  envisioned  through  the  use  of  potential-swing  operation  with  energy  recovery;  however, 
practical  constraints  on  cell  geometry,  aerogel  properties  and  product  concentration  make  this 
specific  application  of  carbon  aerogel  CDI  very  challenging.  The  primary  competing 
technologies  for  BW  desalination  are  electrodialysis  reversal  (EDR)  and  reverse  osmosis  (RO). 
Published  estimates  for  the  desalination  of  1600  ppm  BW  indicates  that  RO  consumes  -8.5  Wh 
gal'1  (8.1  kJ  L'1),  while  ED  consumes  -7.7  Wh  gal'1  (7.3  kJ  L'1).  The  estimated  energy 
requirement  for  desalination  of  2000  ppm  BW  with  carbon  aerogel  CDI  is  -0.53-2.5  Wh  gal’1 
(0.50-2.4  kJ  L’1),  depending  upon  voltage,  flow  rate,  energy  recovery,  electrode  gap,  CAC 
density  and  thickness,  and  other  important  variables.  As  previously  discussed,  the  estimated 
energy  requirements  for  BW  desalination  with  carbon  aerogel  CDI,  RO  and  EDR  are  compared 
in  Table  8. 

Discussion 

Parasitic  Current.  From  the  analysis  presented  in  Tables  6a  through  7i,  it  is  evident  that  the 
open-channel  cell  design  used  in  carbon  aerogel  CDI  minimizes  the  mechanical  energy  required 
for  pumping  to  a  level  where  energy  requirements  are  dominated  by  the  stored  electrical  energy 
and  the  energy  consumed  by  parasitic  electrochemical  reactions  such  as  oxygen  reduction.  Note 
that  the  cell  design  also  minimizes  ohmic  losses  to  a  very  low  level.  Clearly,  efforts  to  improve 
process  efficiency  should  be  directed  towards  the  elimination  or  avoidance  of  parasitic  cell 
current.  It  is  believed  that  this  can  be  done  through  stripping,  catalytic  surfaces  that  eliminate 
reactants  upstream,  the  use  of  suitable  electrode  poisons  (inhibitors),  or  other  conductive  aerogels 
with  greater  activation  overpotentials  for  oxygen  reduction  and  other  relevant  reactions. 
Additional  work  needs  to  be  done  to  quantitatively  determine  parasitic  current  as  a  function  of  all 
relevant  operating  parameters. 

Relative  Challenges  of  BW  and  SW  Applications.  SW  applications  will  require  much  larger 
quantities  of  aerogel  for  the  electrosorption  of  salt  than  BW  applications.  In  this  respect,  SW 
applications  are  expected  to  be  more  difficult.  However,  SW  applications  are  more  tolerant  of 
low  recovery  (product/feed)  ratios,  which  could  be  advantageous.  Recovery  ratios  of  0.50-0.85 
are  required  for  BW  desalination,  whereas  values  as  low  as  0.15  have  been  found  acceptable  for 
some  SW  desalination  processes. 

Scale  Formation.  Most  of  the  salts  found  in  SW  and  BW  have  inverted  solubilities,  which 
means  that  their  solubilities  decrease  with  increasing  temperature.  Thus,  precipitation  and  scale- 
formation  problems  encountered  at  heat  transfer  surfaces  in  thermal  desalination  processes  such 
as  MED-TVC,  MED-MVC  and  MSF  become  worse  as  operating  temperatures  are  increased. 
Furthermore,  the  thermal  decomposition  of  carbonates  such  as  CaC03  produces  C02  and  scale¬ 
forming  hydroxides.  Such  scale  formation  can  be  mitigated  by  suppressing  the  pH  with  HC1  or 
H2S04  additions  to  the  feed.  Unfortunately,  such  acid  dosing  accelerates  corrosion  and  can 
contribute  to  scale  formation  through  the  formation  of  CaS04,  which  is  also  relatively  insoluble. 
By  using  ion  exchange  columns  to  remove  scale-forming  solutes  such  as  Ca2+  and  Mg2+ 
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(hardness)  from  SW  and  BW,  a  relativley  soft  NaCl  solution  is  obtained.  This  softened  stream 
can  then  be  fed  to  the  evaporator  of  a  thermal  desalination  process  with  substantially  less  risk  of 
scale  formation.  This  strategy  has  been  explored  by  Donath  [44].  RO,  ED  and  EDR  are  also 
susceptible  to  scale-formation.  As  previously  discussed,  concentration  polarization  at  the  brine- 
membrane  interface  in  a  RO  system  can  result  in  localized  supersaturation,  precipitation  and 
scale  formation.  Here  too  chemicals  are  added  to  the  feed  to  help  minimize  scale  formation 
problems.  The  ion  separation  inherent  in  carbon  aerogel  CDI  results  in  the  formation  of  high- 
solubility  acids  and  bases  at  anodes  and  cathodes,  respectively,  during  regeneration  (discharge). 
Thus,  scale  formation  due  to  precipitation  from  supersaturated  solution  at  the  electrode¬ 
electrolyte  interface  is  avoided. 

Other  Operational  Advantages  of  Carbon  Aerogel  CDI.  Compressor  maintenance  and  repair 
required  by  MVC  and  MED-MVC  are  avoided.  Since  no  membranes  or  high  pressure  pumps  are 
needed,  carbon  aerogel  CDI  also  offers  operational  advantages  over  RO,  ED  and  EDR.  Standard 
pressure  RO  systems  require  differential  pressures  of  400-650  psi  (27-44  atm),  while  low- 
pressure  systems  require  200-300  psi  (14-20  atm).  The  pressure  drop  across  a  carbon  aerogel 
CDI  system  used  for  BW  desalination  is  expected  to  be  far  less  (Tables  6a  &  6d;  7a  &  7e). 
Typically,  the  shelf  life  of  RO  membranes  is  3  to  5  years,  requiring  cold  storage  during  that 
period  of  time.  Carbon  aerogel  electrodes  have  an  indefinite  shelf  life  and  require  no  cold 
storage.  Furthermore,  RO  membranes  are  susceptible  to  many  of  the  chemicals  used  for  water 
treatment.  Note  that  polyamide  has  been  found  to  be  sensitive  to  attack  by  chlorine,  whereas 
carbon  aerogel  is  believed  to  be  relatively  resistant.  The  cost  of  replacing  damaged  membranes 
is  significant.  Membrane  replacement  accounted  for  -52%  of  the  production  cost  in  early  RO 
systems  [43],  and  for  -32%  in  more  recent  systems  [4],  In  modem  RO  plants,  membrane 
replacement  is  believed  to  account  for  -15%  of  the  production  cost.  These  costs  will  be  driven 
to  higher  values  if  more  expensive  membranes  are  used. 

Recent  Improvements  in  RO  Technology.  Recent  improvements  in  RO  technology  may  lead  to 
improved  energy  efficiency  [47].  The  cost  of  desalinating  1000  gallons  of  oceanic  SW  (-35,000 
ppm)  with  RO  is  given  as  $3.34  by  Shield  and  Moch.  In  their  computational  analysis, 
approximately  40%  of  that  cost  appears  to  be  due  to  the  required  electrical  power.  Since  they 
assume  that  electrical  power  costs  approximately  $0.08  per  kW,  the  corresponding  electrical 
power  requirement  must  have  been  17-18  Wh  gal'1  (16-17  kJ  L'1).  No  estimates  for  BW 
desalination  are  given.  Carbon  aerogel  CDI  could  be  used  as  a  polishing  process  for  low- 
concentration  product  streams  from  such  high-efficiency  processes. 

Selective  Removal  of  Ions.  Ion  exchange  is  now  used  to  remove  heavy  metals  from  contaminated 
ground  water.  Unfortunately,  the  chemical  regeneration  of  saturated  ion  exchange  resin  produces 
a  significant  amount  of  secondary  waste.  A  variety  of  electrochemical  alternatives  to  ion 
exchange  are  now  under  investigation.  For  example,  polypyrolle  films  on  reticulated  vitreous 
carbon  electrodes  are  being  used  for  in  situ  reduction  of  Cr(VI)  [45].  This  approach  appears  to 
suffer  from  a  gradual  degradation  of  the  electrodes  due  to  a  loss  of  polypyrolle.  Another  process 
involves  the  use  of  electrodes  coated  with  films  of  electroactive  ferrocyanide  for  selective 
removal  Cs  from  solutions  of  sodium  salts  [46].  Here  too  electrode  life  may  be  limited  by  the 
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stability  of  the  electroactive  film.  Carbon  aerogel  electrodes  appear  to  have  the  necessary 
selectivity  and  stability  to  enable  them  to  serve  as  a  viable  alternative  to  ion  exchange. 
Treatability  tests  at  LLNL  have  shown  that  35  ppm  Cr(VI)  can  be  selectively  removed  from 
contaminated  ground  water  with  a  530  ppm  TDS  background  [24.27],  Reversibility  was 
demonstrated  by  quantitiative  recovery  of  all  electrosorbed  chromium  during  regeneration. 
Carbon  aerogel  CDI  has  also  been  used  to  remove  Cu,  Zn,  Cd,  Pb,  Cr,  Mn,  Co  and  U  from 
process  solutions  and  natural  waters  by  electrosorption  on  CAC  electrodes  [24].  This  type  of 
application  was  the  incentive  for  development  of  carbon  aerogel  CDI. 

Summary 

The  minimum  theoretical  work  required  by  an  isothermal  process  to  separate  sea  water  (-35,000 
ppm  NaCl)  or  brackish  water  (-1000-10,000  ppm  NaCl)  into  product  and  concentrate  streams 
has  been  calculated.  This  calculation  requires  that  the  activity  coefficient  of  NaCl  be  estimated 
with  the  Debye-Huckel  Theory.  Assuming  a  recovery  ratio  of  0.5,  the  minimum  theoretical  work 
required  for  separation  of  SW  (35,000  ppm  NaCl)  into  product  (100  ppm  NaCl)  and  concentrate 
streams  is  -2.0  Wh  gal’1  (1.9  kJ  L’1),  whereas  the  minimum  theoretical  work  for  generating  the 
same  product  from  BW  (2000  ppm  NaCl)  is  only  -0.10  Wh  gal 1  (0.09  kJ  L  ').  Data  from 
several  operating  desalination  processes  were  published  by  Wade  in  1993  [2].  His  work 
indicates  that  RO  consumes  -25-36  Wh  gal’1  (24-34  kJ  L'1),  MVC  consumes  -30-41  Wh  gal 1 
(29-39  kJ  L'1),  ED  consumes  -60-75  Wh  gal’1  (57-71  kJ  L'1),  MED-TVC  consumes  -215-315 
Wh  gal’1  (205-300  kJ  L'1),  and  MSF  consumes  -320  Wh  gal’1  (304  kJ  L’1).  Simple  evaporation 
without  energy  recovery  consumes  -2720  Wh  gal’1  (2587  kJ  L’1).  More  recent  computational 
analysis  by  Shields  and  Moch  assumes  that  RO  requires  17-18  Wh  gal’1  (16-17  kJ  L’1)  [47].  The 
intrinsic  energy  requirement  for  SW  desalination  with  carbon  aerogel  CDI  is  QV/2,  which  is  -36- 
37  Wh  gal’1  (34-35  kJ  L’1).  An  energy  requirement  less  than  that  needed  for  RO  can  be 
envisioned  through  the  use  of  potential-swing  operation  with  energy  recovery;  however,  practical 
constraints  on  cell  geometry,  aerogel  properties  and  product  concentration  make  this  specific 
application  of  carbon  aerogel  CDI  very  challenging.  The  primary  competing  technologies  for 
BW  desalination  are  electrodialysis  reversal  (EDR)  and  reverse  osmosis  (RO)  [3a,3b,4]. 
Published  data  for  the  desalination  of  1600  ppm  BW  indicates  that  RO  consumes  -8.5  Wh  gal’1 
(8.1  kJ  L’1),  while  ED  consumes  -7.7  Wh  gal'1  (7.3  kJ  L’1).  The  estimated  energy  requirement 
for  desalination  of  2000  ppm  BW  with  carbon  aerogel  CDI  is  -0.53-2.5  Wh  gal’1  (0.50-2.4  kJ 
L’1),  depending  upon  voltage,  flow  rate,  energy  recovery,  electrode  gap,  CAC  density  and 
thickness,  and  other  important  variables.  Attractive  scenarios  for  BW  desalination  have  been 
found;  however,  cell  designs  have  demanding  tolerances.  With  additional  investments  in 
research  and  development,  it  may  be  possible  to  develop  carbon  aerogel  CDI  to  a  point  where 
practical,  large-scale  desalination  of  BW  is  possible.  However,  such  applications  will  probably 
require  significant  reductions  in  the  costs  of  the  aerogel  materials. 

Carbon  aerogel  electrodes  have  demonstrated  the  necessary  selectivity  and  stability  to  enable 
them  to  serve  as  a  viable  alternative  to  ion  exchange  resins  for  the  remediation  of  Cr(VI)- 
contaminated  ground  water.  Field  tests  at  LLNL  have  shown  that  carbon  aerogel  electrodes  can 
be  used  to  selectively  remove  35  ppb  Cr(VI)  from  contaminated  ground  water  with  a  530  ppm 
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TDS  background.  Consequently,  it  may  be  possible  to  avoid  the  generation  of  large  quantities  of 
secondary  waste,  produced  during  chemical  regeneration  of  ion  exchange  columns. 

Future  Work 

Since  numerous  assumptions  were  made  to  render  computations  tractable,  the  estimates 
presented  here  are  imperfect.  As  knowledge  of  this  novel  aerogel-based  process  improves, 
calculations  should  be  refined.  Furthermore,  experimental  data  from  experiments  with 
concentrated  NaCl  and  actual  BW  should  be  used  to  verify  predictions. 

Estimation  of  the  minimum  product  concentration,  CP ,  should  be  improved.  These  calculations 
evaluate  the  area-specific  electrosorption  capacity,  cr,  at  the  feed  concentration,  CF.  Since  the 
electrodes  become  saturated  with  salt  at  the  point  when  equilibrium  is  established  between  the 
feed  stream  and  the  carbon  aerogel,  this  seems  like  a  reasonable  assumption.  However, 
immediately  following  polarization,  the  equilibrium  established  at  the  stack  outlet  is  at  Cf  instead 
of  CF.  In  the  future,  any  equation  used  to  estimate  the  minimum  product  concentration  should 
use  a  value  of  cr  evaluated  at  the  outlet  conditions. 

Ultimately,  the  concentration  profile  in  the  stack  of  carbon  aerogel  electrodes  should  be 
calculated  by  analytically  or  numerically  integrating  the  first-order  hyperbolic  partical 
differential  equations  that  describe  the  system.  The  analytical  methods  for  chromatographic 
systems  outlined  by  Rhee,  Aris  and  Amundsen  are  believed  to  be  applicable  [48],  Chiba  has 
solved  the  equations  for  a  carbon  aerogel  CDI  system  numerically  [49].  Both  sources  should  be 
consulted. 

To  overcome  limitations  in  product  concentration,  parametric  pumping  should  be  explored.  For 
example,  feed  should  be  admitted  after  polarization  so  that  a  steep  concentration  gradient  can  be 
established  along  the  length  of  the  channel  (stratification  of  the  column).  Subsequently,  the 
electrodes  should  be  electrically  discharged  (regenerated)  with  zero  flow  (stagnant).  After 
discharge,  the  concentrate  should  be  removed  from  the  high-concentration  end  of  the  stack, 
where  the  feed  is  introduced.  Such  reverse  flow  (backflushing)  should  enable  the  steep 
concentration  gradient  to  be  maintained.  Gas  displacement  may  be  required.  Some  aspects  of 
this  mode  of  operation  have  already  been  explored  by  Oren  and  Soffer  [13-15]. 
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Table  1.  Estimates  of  Total  Energy  Requirements  for  SW  Desalination 


Note:  MSF  (multistage  flash);  MED  (multiple  effect  distillation);  TVC  (thermo  vapor  compression):  MVC 
(mechanical  vapor  compression);  RO  (reverse  osmosis);  ERT  (energy  recovery  turbines). 


Source:  N.  M.  Wade,  "Technical  and  Economic  Evaluation  of  Distillation  and  Reverse  Osmosis  Desalination 
Processes,"  Desalination,  Vol.  93,  pp.  343-363,  (1993). 

Table  2.  Heuristic  Estimates  of  Energy  Requirements  for  Reverse  Osmosis 


Salt 

Concentration 

Osmotic 

Pressure 

System 

Category 

Pressure 

(Low) 

Pressure 

(High) 

Energy 

Requirement 

(Low) 

Energy  Requirement 
(High) 

PPtn 

psi 

psi 

psi 

Wh  gal" 

Wh  gaf  * 

1000 

10 

LP 

200 

300 

2.06 

3.09 

1500 

15 

LP 

200 

300 

2.06 

3.09 

2000 

20 

LP 

200 

300 

2.06 

3.09 

2500 

25 

LP 

200 

300 

2.06 

3.09 

3000 

30 

LP 

200 

300 

2.06 

3.09 

3500 

35 

LP/SP 

200 

650 

2.06 

6.69 

4000 

40 

SP 

400 

650 

4.11 

6.69 

4500 

45 

SP 

400 

650 

4.11 

6.69 

5000 

50 

SP 

400 

650 

4.11 

6.69 

10000 

100 

SP 

400 

650 

4.1 1 

6.69 

35000 

350 

SW 

800 

1500 

8.23 

15.43 

Note:  LP  (low  pressure);  SP  (standard  pressure);  SW  (sea  water).  The  energy  requirement  is  calculated  from  the 
pressure  alone  and  is  based  on  a  unit  volume  of  product,  with  an  assumed  recovery  ratio  of  0. 70. 

Sources:  W.  J.  Conlon,  "Membrane  Processes,"  Chapt.  II,  Water  Quality  and  Treatment.  A  Handbook  of 
Community  Water  Supplies.  4th  Ed.,  R.  W.  Pontius,  Ed.,  American  Water  Works  Association,  McGraw-Hill,  San 
Francisco,  CA,  pp.  709-746  (1990);  R.  M.  Clark,  "Water  Supply,  "  Chapt.  5.  Standard  Handbook  of  Environmental 
Engineering.  R.  A.  Corbitt,  Ed.,  McGraw-Hill,  San  Francisco,  CA,  pp.  5.146-149  (1989). 
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Table  3.  Estimates  of  Energy  Requirements  for  Electrodialysis  Reversal  (EDR) 


Salt 

Concentration 

EDR  Membrane  Stack 

EDR  Pumping 

EDR  Total 

ppm 

Wh  gaf 

igigEESlBil 

1000 

2 

2.5 

4.73 

1500 

3 

2.5 

5.78 

2000 

4 

2.5 

6.83 

2500 

5 

2.5 

7.88 

3000 

6 

2.5 

8.93 

3500 

7 

2.5 

9.98 

4000 

8 

2.5 

11.03 

4500 

9 

2.5 

12.08 

5000 

10 

2.5 

13.13 

10000 

20 

2.5 

23.63 

35000 

70 

2.5 

76.13 

Source:  W.  J.  Conlon,  "Membrane  Processes,  ”  Chapt.  I  I,  Water  Oualit\’  and  Treatment.  A  Handbook  of 
Community  Water  Supplies.  4th  Ed.,  R.  W.  Pontius,  Ed.,  American  Water  Works  Association,  McGraw-Hill,  San 
Francisco,  CA,  pp.  737  (1990).  The  total  assumes  5%  is  due  to  instrumentation. 


Table  4.  Estimated  Energy  Requirement  for  City  of  Melville  Desalting  Plant 


Plant  Feature 

RO 

EDR 

Well  Concentration  (TDS,  ppm) 

1600 

1600 

Feed  Concentration  (TDS,  ppm) 

1630 

1600 

Concentrate  Concentration  (TDS,  ppm) 

5300 

7500 

Product  Concentration  (TDS,  ppm) 

60 

490 

Product-Feed  (Recovery)  Ratio 

0.700 

0.842 

Feed  Capacity  (L  s*1) 

23.9 

26.5 

Product  Capacity  (L  s'1) 

16.7 

22.3 

Supply  Capacity  (L  s'1) 

22.3 

22.3 

Blending  Ratio  (Product:  Well) 

3:1 

none 

Product  Concentration  After  Blending  (TDS,  ppm) 

445 

490 

Operating  Pressure  (kPa/psi) 

1,860/270 

380/55 

Energy  Requirement  for  Pressure  Drop  Alone  (Wh  gal1) 

2.78 

0.47 

Capital  Cost  ($K  Canadian) 

1,965 

1,942 

Operating  Costs  ($  mJ  product) 

Electricity  ($0,049  kWh) 

0.11 

0.10 

Chemicals 

0.12 

0.02 

Membrane  Replacement 

0.08 

0.06 

Operation  &  Maintenance 

0.04 

0.03 

Filters,  Materials  &  Spare  Parts 

0.03 

0.03 

Total  Operating  Cost  -No  Blending  ($  mJ  product) 

0.38 

0.24 

Total  Operating  Cost  -  With  Blending  ($  mJ  product) 

0.29 

0.24 

Energy  Requirement  Based  on  Electricity  (Wh  gaf1) 

8.49 

7.71 

Source:  R.  C.  Harries,  D.  Elyanow,  D.  N.  Heshka,  K.  L.  Fischer,  “Desalination  of  Brackish  Groundwater  for  a 
Prarie  Community  Using  Electrodialysis  Reversal,  "  Desalination,  Vol.  84.  pp.  109-121,  (1991).  The  value  of  7. 71 
Wh  gal1  for  EDR  is  higher  than  that  estimated  with  the  heuristic  rule,  which  is  5.99  Wh  gaf' . 
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Table  6b.  Electrode  Gap  (in)  vs.  Feed  Concentration  for  Conditions  in  Table  6a. 


Table  6c.  Channel  Length  (ft)  vs.  Feed  Concentration  for  Conditions  in  Table  6a. 


Table  6d.  Pressure  Drop  (psi)  vs.  Feed  Concentration  for  Conditions  in  Table  6a. 
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Table  6e.  Energy  Required  (Wh  gal'1)  vs.  Feed  Concentration  -  Table  6a  -  No  Depletion. 


Feed 

Cone. 

(ppm) 

Case  1: 
1.2V&  100 
ml  min’1  stack’1 

Case  2: 

1.2  V  &  25 
ml  min'1  stack'1 

Case  3: 
1.2V&  10 
ml  min'1  stack'1 

Case  6: 
increased 
CAC  density 

Stored  Energy 
-1/2  QV 
@  1.2  V 

Minimum 

Theoretical 

@R=0.5 

1000 

2.211 

1.657 

1.395 

0.538 

0.937 

0.0412 

1500 

3.165 

2.377 

0.640 

1.458 

0.0677 

2000 

4.218 

3.054 

2.592 

0.770 

1.979 

0.0946 

2500 

11.022 

3.701 

3.152 

0.914 

0.1217 

3000 

4.329 

3.694 

3.021 

0.1488 

3500 

lillillitiilll 

4.949 

4.224 

3.541 

0.1760 

4000 

5.642 

4.743 

4.062 

0.2032 

4500 

10.580 

5.253 

4.583 

0.2305 

5000 

5.757 

5.104 

0.2579 

10000 

10.3n 

0.5338 

35000 

IBfBllIlB 

36350 

1.9895 

Table  6f.  Energy  Required  (Wh  gal'1)  vs.  Feed  Concentration  -  Table  6a  -  With  Depletion. 


Feed 

Cone. 

(ppm) 

Case  1 : 
1.2V&  100 
ml  min'1  stack'1 

Case  2: 

1.2  V  &  25 
ml  min'1  stack'1 

Case  3: 
1.2V&  10 
ml  min'1  stack'1 

Case  6: 
increased 
CAC  density 

Stored  Energy 
-1/2  QV 
@1.2  V 

Minimum 
Theoretical 
@  R=0.5 

1000 

0.693 

1.027 

1.445 

0.543 

0.937 

0.0412 

1500 

0.919 

1.350 

2.099 

0.707 

1.458 

0.0677 

2000 

1.343 

1.613 

2.510 

1.01 1 

1.979 

0.0946 

2500 

7.577 

1.873 

2.802 

1.806 

2.500 

0.1217 

3000 

2.135 

3.065 

m— 

3.021 

0.1488 

3500 

2.408 

3.325 

3.541 

0.1760 

4000 

2.766 

3.585 

4.062 

0.2032 

4500 

7.382 

3.845 

4.583 

0.2305 

5000 

4.107 

5.104 

0.2579 

10000 

10.311 

0.5338 

35000 

36.350 

1.9895 
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found.  Computational  results  for  other  feed  concentrations  are  summarized  in  Tables  7  b  through  7i,  as  well  as  in  Figures  7  a  through  J. 
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Table  7b.  Minimum  Product  Concentration  (ppm)  vs.  Feed  Concentration  -  Table  7a. 


Feed 

Cone. 

(ppm) 

Case  1: 

1.2  V&  100 
ml  min"1  stack'1 

Case  2: 

1.2  V  &  25 
ml  min'1  stack'1 

Case  3: 
1.2V&  10 
ml  min'1  stack'1 

Case  4: 

0.6  V&  10 
ml  min'1  stack'1 

Case  5: 

Case  3  with 
smaller  gap 

Case  6: 
increased 
CAC  density 

1000 

211 

226 

194 

770 

182 

173 

1500 

618 

203 

211 

1243 

179 

186 

2000 

1045 

596 

189 

1722 

180 

174 

2500 

1485 

1007 

574 

2204 

185 

171 

3000 

1933 

1431 

975 

2689 

172 

178 

3500 

2387 

1863 

1387 

3175 

199 

167 

4000 

2845 

2301 

1808 

3663 

407 

189 

4500 

3307 

2746 

2236 

4152 

789 

182 

5000 

3772 

3194 

2670 

4642 

1180 

199 

10000 

8514 

7815 

7181 

9567 

5378 

191 

35000 

32903 

31917 

31022 

34389 

28478 

9742 

Table  7c.  Maximum  Recovery  (Product/Feed)  Ratio  vs.  Feed  Concentration  -  Table  7a. 


Feed 

Case  1: 

Case  2: 

Case  3: 

Case  4: 

Case  5: 

Case  6: 

Cone. 

1.2  V&  100 

1.2  V  &  25 

1.2  V  &  10 

0.6  V&  10 

Case  3  with 

increased 

(ppm) 

ml  min'1  stack'1 

ml  min'1  stack'1 

ml  min'1  stack'1 

ml  min'1  stack  1 

smaller  gap 

CAC  density 

1000 

0.500 

0.600 

0.650 

0.500 

0.750 

0.920 

1500 

0.500 

0.500 

0.565 

0.500 

0.675 

2000 

0.500 

0.500 

0.500 

0.500 

0.620 

2500 

0.500 

0.500 

0.500 

0.500 

0.577 

0.840 

3000 

0.500 

0.500 

0.500 

0.500 

0.540 

0.820 

3500 

0.500 

0.500 

0.500 

0.500 

0.512 

0.801 

4000 

0.500 

0.500 

0.500 

0.500 

0.500 

0.785 

4500 

0.500 

0.500 

0.500 

0.500 

0.500 

0.769 

5000 

0.500 

0.500 

0.500 

0.500 

0.500 

0.755 

10000 

0.500 

0.500 

0.500 

0.500 

0.500 

0.646 

35000 

0.500 

0.500 

0.500 

0.500 

0.500 

0.500 
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Table  7d.  Channel  Length  (ft)  vs.  Feed  Concentration  for  Conditions  in  Table  7a. 


Feed 

Cone. 

(ppm) 

Case  1 : 

1.2  V  &  100 
ml  min’1  stack'1 

Case  2: 

1.2  V  &  25 
ml  min'1  stack’1 

Case  3: 

1.2  V&  10 
ml  min'1  stack'1 

Case  4: 
0.6V&  10 
ml  min’1  stack’1 

Case  5: 
Case  3  with 
smaller  gap 

Case  6: 
increased 
CAC  density 

1000 

346.654 

57.776 

18.666 

34.665 

18.944 

0.797 

1500 

346.654 

86.663 

26.689 

34.665 

27.364 

1.133 

2000 

346.654 

86.663 

34.665 

34.665 

34.833 

1.455 

2500 

346.654 

86.663 

34.665 

34.665 

41.664 

1.745 

3000 

346.654 

86.663 

34.665 

34.665 

48.413 

2.011 

3500 

346.654 

86.663 

34.665 

34.665 

54.168 

2.277 

4000 

346.654 

86.663 

34.665 

34.665 

56.832 

2.510 

4500 

346.654 

86.663 

34.665 

34.665 

56.832 

2.753 

5000 

346.654 

86.663 

34.665 

34.665 

56.832 

2.974 

10000 

346.654 

86.663 

34.665 

34.665 

56.832 

5.021 

35000 

346.654 

86.663 

34.665 

34.665 

56.832 

9.163 

Table  7e.  Pressure  Drop  (psi)  vs.  Feed  Concentration  for  Conditions  in  Table  7a. 


Feed 

Cone. 

(ppm) 

Case  1: 

1.2  V&  100 
ml  min’1  stack’1 

Case  2: 

1.2  V  &  25 
ml  min’1  stack  1 

Case  3: 

1.2  V&  10 
ml  min'1  stack’1 

Case  4: 

0.6V&  10 
ml  min’1  stack’1 

Case  5: 

Case  3  with 
smaller  gap 

Case  6: 

increased 
CAC  density 

1000 

3.52 

fgBKHHI 

iSEMEH 

1500 

3.52 

2000 

3.52 

■BW 

wmmm 

KgfESEgjH 

2500 

3.52 

WKE5SSMM 

1.89x1  O^ 

■bhh 

KSBEHH 

3000 

3.52 

■iggiEgH 

3500 

3.52 

W/EE5MM 

HEBiM 

1.89x10 

4000 

3.52 

WHESSMM 

1. 89x10’'“ 

1 .89x1 0'^ 

UfeSLlH 

isiiiHi 

4500 

3.52 

1.89x10^ 

msmmm 

i— 

5000 

3.52 

mmbbemm 

1 .89x10’'“ 

;1WI 

urmur 

msesehm 

10000 

3.52 

1.89x10^ 

1 .89x10’'“ 

\W3E3MSi 

35000 

3.52 

iheqeim 

1.89x10'^ 

1.89x10^ 

\wmvm 

38 


Lawrence  Livermore  National  Laboratory  Report  No.  UCRL-lD-l 25298  Rev.  1 


Table  7f.  Minimum  Energy  Requirement  (Wh  gal'1)  vs.  Feed  Concentration  -  Table  7a. 


Feed 

Cone. 

(ppm) 

Case  l: 

1.2V&  100 
ml  min  '  stack'1 

Case  2: 

1.2  V  &  25 
ml  min'1  stack'1 

Case  3: 
1.2V&  10 
ml  min'1  stack'1 

Case  4: 

0.6  V&  10 
ml  min’1  stack’1 

Case  5: 
Case  3  with 
smaller  gap 

Case  6: 
increased 
CAC  density 

1000 

0.030 

0.032 

0.038 

0.002 

0.046 

0.076 

1500 

0.055 

0.002 

0.072 

0.112 

2000 

0.019 

0.044 

0.002 

0.096 

0.149 

2500 

0.017 

0.039 

0.068 

0.001 

0.118 

0.183 

3000 

0.016 

0.035 

0.061 

0.001 

0.143 

0.214 

3500 

0.015 

0.032 

0.055 

0.001 

0.162 

0.247 

4000 

0.014 

0.030 

0.051 

0.001 

0.159 

0.274 

4500 

0.013 

0.028 

0.048 

0.001 

0.144 

0.305 

5000 

0.012 

0.027 

0.045 

0.001 

0.133 

0.331 

10000 

0.019 

0.032 

0.001 

0.089 

0.603 

35000 

0.055 

0.012 

0.020 

0.000 

0.053 

0.859 

Table  7g.  Stored  Electrical  Energy  ~  1/2  QV  (Wh  gal'1)  vs.  Feed  Concentration  -  Table  7a. 


Feed 

Case  1 : 

Case  2: 

Case  3: 

Case  4: 

Case  5: 

Case  6: 

Cone. 

1.2V&  100 

1.2  V&25 

1.2V&  10 

0.6V&  10 

Case  3  with 

increased 

(ppm) 

ml  min'1  stack'1 

ml  min'1  stack'1 

ml  min'1  stack'1 

ml  min'1  stack  1 

smaller  gap 

CAC  density 

1000 

0.822 

0.806 

0.840 

0.120 

0.852 

0.9328 

1500 

0.919 

1.351 

1.342 

0.134 

1.376 

1.4821 

2000 

0.994 

1.462 

1.887 

0.145 

1.896 

2.0603 

2500 

1.057 

1.555 

.2.006 

0.154 

2.411 

2.6284 

3000 

1.112 

1.635 

2.109 

0.162 

2.946 

3.1847 

3500 

1.160 

1.705 

2.200 

0.169 

3.438 

3.7604 

4000 

1.203 

1.769 

2.283 

0.175 

3.742 

4.3005 

4500 

1.243 

1.827 

2.358 

0.181 

3.865 

4.8718 

5000 

1.279 

1.881 

2.427 

0.186 

3.979 

5.4175 

10000 

1.548 

2.276 

2.936 

0.226 

4.814 

11.0685 

35000 

2.184 

3.211 

4.143 

0.318 

6.793 

28.5003 
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Table  7h.  Energy  Required  (Wh  gal'1)  vs.  Feed  Concentration  -  Table  7a  -  No  Depletion. 


Feed 

Case  1: 

Case  2: 

Case  3: 

Case  4: 

Case  5: 

Case  6: 

Cone. 

1.2  V&  100 

1.2  V  &  25 

1.2V&  10 

0.6  V&  10 

Case  3  with 

increased 

(ppm) 

ml  min’1  stack*1 

ml  min"'  stack"1 

ml  min"1  stack*1 

ml  min*1  stack'1 

smaller  gap 

CAC  density 

1000 

1.939 

1.422 

1.245 

0.369 

1.260 

0.3461 

1500 

1.986 

2.203 

1.849 

0.376 

1.892 

0.5320 

2000 

2.023 

2.258 

2.472 

0.381 

2.480 

0.7249 

2500 

2.054 

2.304 

2.531 

0.386 

3.038 

0.9123 

3000 

2.081 

2.343 

2.582 

0.389 

3.602 

1.0943 

3500 

2.105 

2.378 

2.627 

0.393 

4.101 

1.2820 

4000 

2.127 

2.410 

2.668 

0.396 

4.370 

1.4568 

4500 

2.146 

2.439 

2.705 

0.399 

4.431 

1.6416 

5000 

2.164 

2.466 

2.740 

0.402 

4.488 

1.8173 

10000 

2.298 

2.662 

2.993 

0.421 

4.905 

3.6217 

35000 

2.616 

3.129 

3.595 

0.467 

5.893 

9.0503 

Table  7i.  Energy  Required  (Wh  gal'1)  vs.  Feed  Concentration  -  Table  7a  -  With  Depletion. 


Feed 

Cone. 

(ppm) 

Case  1: 
1.2V&  100 
ml  min*1  stack'1 

Case  2: 

1.2  V&25 
ml  min*1  stack'1 

Case  3: 
1.2V&  10 
ml  min'1  stack'1 

Case  4: 

0.6  V&  10 
ml  min'1  stack  1 

Case  5: 
Case  3  with 
smaller  gap 

Case  6: 
increased 
CAC  density 

1000 

0.657 

0.970 

1.716 

0.809 

1.969 

0.5477 

1500 

0.705 

1.290 

2.107 

0.815 

2.241 

0.8060 

0.742 

1.345 

2.444 

0.821 

2.502 

1.0634 

0.774 

1.391 

2.503 

0.825 

2.760 

1.3048 

3000 

0.801 

1.431 

2.554 

0.829 

3.027 

1.5330 

3500 

0.825 

1.466 

2.599 

0.832 

3.273 

1.7639 

4000 

0.847 

1.498 

2.640 

0.836 

3.425 

1.9744 

4500 

0.866 

1.527 

2.677 

0.838 

3.487 

2.1943 

5000 

0.884 

1.554 

2.712 

0.841 

3.543 

2.4000 

10000 

1.019 

1.750 

2.965 

0.860 

3.960 

4.4130 

35000 

1.336 

2.218 

3.568 

0.907 

4.949 

9.9963 

Table  8.  A  Comparison  of  Desalination  Processes  for  Brackish  Water 


EDR 

RO 

CAC  CD1 

CAC  CD1 

Source 

Tables  3  &  4 

Tables  2  &  4 

Table  7  -  Case  3 

Table  7  -  Case  6 

Feed  Concentration  (ppm) 

1600 

1630 

1500/2000 

1500/2000 

Product  Concentration  (ppm) 

490 

60 

211/189 

186/174 

Recovery  (Product/Feed)  Ratio 

0.842 

0.700 

0.565/0.500 

0.890/0.863 

Pressure  Differential  (psi) 

55 

270 

■lymnrrmM 

wjvmwm 

Energy  Requirement  (Wh  gal'1)  -  Lowest 

5.78-6.83 

2.06-3.09 

1.85/2.47 

0.53/0.73 

Energy  Requirement  (Wh  gal'1)  -  Highest 

7.71 

8.49 

2.11/2.44 

0.81/1.06 
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Figure  Captions 

Figure  1.  Electrosorption  of  Na+ and  Cl  ions  on  54  double-sided  CAC  electrodes.  Results  for 
100,  500  and  1000  mS  cm'1. 

Figure  2.  Electrosorption  of  Na+  and  Cl-  on  100  double-sided  electrodes.  Results  for  0.6  and 
1.2  V  compared. 

Figure  3.  Electrosorption  capacity  for  NaCl  on  carbon  aerogel  composite  (CAC)  electrodes. 
Comparison  of  first  regression  analysis  (Equation  5)  to  experimental  data. 

Figure  4.  Electrosorption  capacity  for  NaCl  on  carbon  aerogel  composite  (CAC)  electrodes. 
Comparison  of  second  regression  analysis  (Equation  6)  to  experimental  data. 

Figure  5.  Moody  friction  factor,  based  upon  measured  pressure  drop,  correlated  with  Reynolds 
number. 

Figure6a.  Assumed  recovery  ratios  (0.50  or  0.85)  at  various  feed  concentrations,  for  Cases  1 
through  6  given  in  Table  6a. 

Figure  6b.  Calculated  electrode  gap  (channel  height)  at  various  feed  concentrations,  for  Cases  1 
through  6  given  in  Table  6a.  Note  that  Table  6a  only  provides  the  value  at  2000  ppm. 

Figure  6c.  Calculated  energy  requirements  at  various  feed  concentrations,  for  Cases  1  through  6 
given  in  Table  6a.  These  estimates  assume  no  depletion  of  electroactive  species 
responsible  for  parasitic  current  and  are  therefore  considered  a  “worst  case”  scenario. 

Figure  6d.  Calculated  energy  requirements  shown  in  Figure  6c  with  magnified  scale. 

Figure  6e.  Calculated  energy  reuirements  at  various  feed  concentrations,  for  Cases  1  through  6 
given  in  Table  6a.  These  estimates  assume  that  electroactive  species  responsible  for 
parasitic  current  are  introduced  in  the  feed  and  are  therefore  depleted,  a  “best  case” 
scenario. 

Figure  6f.  Calculated  energy  requirements  shown  in  Figure  6e  with  magnified  scale. 

Figure  7a.  Estimates  of  the  maximum  possible  recovery  (product/feed)  ratios  at  various  feed 
concentrations,  for  cases  1  through  6  given  in  Table  7a. 

Figure  7b.  Estimates  of  the  minimum  possible  salt  concentration  in  the  product  at  various  feed 
concentrations,  for  Cases  1  through  6  given  in  Table  7a. 

Figure  7c.  Calculated  theoretical  minimum  energy  requirement  for  separation  of  feed  into 

product  and  concentrate.  Values  at  various  feed  concentrations,  for  Cases  1  through  6 
given  in  Table  7a. 
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Figure  7d.  Calculated  energy  storage  during  separation  of  feed  into  product  and  concentrate. 
Values  at  various  feed  concentrations,  for  Cases  1  through  6  given  in  Table  7a. 

Figure  7e.  Calculated  energy  requirements  at  various  feed  concentrations,  for  cases  1  through  6 
given  in  Table  7a.  These  estimates  assume  that  electroactive  species  responsible  for 
parasitic  current  are  introduced  in  the  feed  and  are  therefore  depleted,  a  “best  case” 
scenario. 

Figure  If.  Calculated  energy  requirement  at  various  feed  concentrations,  for  Cases  1  through  6 
given  in  Table  6a.  These  estimates  assume  no  depletion  of  electroactive  species 
responsible  for  parasitic  current  and  are  therefore  considered  a  “worst  case”  scenario. 
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